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''Thia  report  discusses  the  possibility  of  extracting  sleep  information  from 
heart  rate  data-  Eight  hours  of  sleep  EEG,  EOG  and  electrocardiograms  were  recorded 
on  FM  magnetic  tape  for  two  nights  for  each  of  eight  subjects.  The  time  in  milli¬ 
seconds  between  heart  beats  was  written  on  digital  magnetic  tape.  Tne  data  were 
grouped  into  records  containing  128  consecutive  beat-to-bcat  intervals  and  a  large 
number  of  descriptors  was  computed  for  each  record.  These  descriptors  for  each 
record  were  the  mean  value,  the  sample  variance,  the  nine- interval  histogram  of  2 
scores,  and  a  set  of  Fourier  transforms.  Analysis  of  variance  was  used  as  a  general 
guide  to  descriptor  significance  for  each  subject.  The  discriminant  analysis  pro¬ 
cedure  described  by  Rao  and  popularized  by  Cooley  and  Lohnes  was  used  to  sleep  stage 
classify  heart  rate  data.  Accuracy  of  the  procedure  was  determined  in  terms  of 
percent  correct  classifications,  correlation  coefficient  of  the  computerized  sleep 
pattern  with  respect  to  the  EEG  hand  scored  pattern,  and  an  empirically  derived  cost 
function. 

The  results  of  this  study  suggest  that  for  a  single  night  of  sleep  a  reasonable 
accuracy  of  sleep  ctagz  classification  is  possiole.  However  the  variability  in 
heart  rate  frerj  night-to-night  for  any  one  individual  produces  unacceptably  poor 
classification  results  on  the  second  night. 
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ABSTRACT 


This  report  discusses  the  possibility  of  extracting  sleep  information 
from  heart  rate  data.  The  recognition  of  sleep  stages  or  even  the  ability  to 
differentiate  sleep  from  wakefulness  using  heart  rate  information  alone  rather 
than  the  conventional  EEG  measures  could  expend  the  scope  of  sleep  studies. 

In  many  situations  where  it  is  desirable  to  evaluate  wakefulness  and  sleep, 

EEG  electrodes  become  unreliable  after  a  few  days  and  the  time  bandwidth 
requirements  of  recording  and  transmitting  the  EEG  are  excessive. 

Eight  hours  cf  sleep  EEG,  EOG  and  electrocardiograms  were  recorded  on 
FM  magnetic  tape  for  two  nights.  The  method  of  data  collection  and  sleep 
scoring  cf  the  EEG  was  reported  in  detail  by  Lessard,  Ford  and  Hughes  erf 
the  USAF  School  of  Aerospace  Medicine  (14).  Copies  of  the  FM  magnetic  tapes 
containing  sleep  data  for  ten  subjects  were  supplied  to  The  University  of  Texas 
Bio-Medical  Engineering  Laboratory  by  the  U.  S.  Air  Force  School  of  Aerospace 
Medicine. 

Double  differentiation  of  the  filtered  electrocardiogram  and  threshold  logic 

units  were  used  to  detect  the  peak  cf  each  R-wave.  The  time  in  milliseconds 

between  heart  beats  was  written  on  digital  magnetic  tape.  The  data  were 

grouped  into  records  containing  128  consecutive  beat-to-beat  intervals  and 

eleven  descriptors  were  computed  for  each  record .  These  descriptors  for  each 

—  2 

record  were  the  mean  value  X,  the  sample  variance  S  ,  and  the  nine-interval 
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histogram  of  the  beat-to-beat  R-R  intervals.  We  represented  each  R-R  inter¬ 
val  as  and  the  mean  value  and  standard  deviation  of  the  1 28  R-R  intervals 
as  X  and  S,  respectively;  then  the  standardised  valu  lor  the  ith heart  beat  was 


The  distribution  of  Z  scores  had  a  mean  value  of  zero  and  standard  deviation 
of  one.  Histogram  intervals  were  one  standard  deviation  wide  and  the  number 
of  standardized  scores  falling  in  each  of  the  nine  intervals  were  used  as  the 
value  of  the  histogram  descriptors.  Each  group  of  128  heart  intervals  hod  a 
different  mean  value,  X  and  sample  standard  deviation,  S. 

The  outer  Intervals  of  the  histogram  extended  from  -•  to  -1-3/4  and  +  1-3/4 
to  +  •. 

Analysis  of  variance  was  used  to  determine  descriptor  significance  for 
each  subject.  This  procedure  tested  the  null  hypothesis  for  each  descriptor, 
fn  other  words,  what  was  the  probability  that  a  descriptor  mean  value  for  awake 
and  the  five  stages  of  sleep  were  equal?  Our  program  used  an  F-tesi  to  com¬ 
pute  the  probability  that  the  aiean  values  of  the  descriptors  for  each  stage  of 
sleep  were  equal.  The  hypothesis  was  rejected  if  P  <  .01 . 

The  discriminant  analysis  procedure  described  by  Rac  (1 7)  and  popularized 
by  Cooley  and  Lohnes  \7)  was  used  to  sleep  stage  classify  heart  rate  data. 
Approximately  one-half  of  the  first  recorded  night  of  sleep  for  each  individual 


lii 


was  used  as  a  training  set  in  the  discriminant  analysis  procedures.  Once  the 
training  set  had  been  obtained  both  nights  of  sleep  for  the  individual  were 
sleep  stage  classified  into  awake  and  stage  1,  2,  3,  4,  and  REM  sleep. 

Accuracy  of  the  procedure  was  determined  in  terms  of  percent  correct  classifi¬ 
cations  ,  correlation  coefficient  of  the  computerized  sleep  pattern  with  respect 
to  the  EEG  hand  scored  pattern  and  an  empirically  derived  cost  function. 

The  results  of  this  study  suggest  that  for  a  single  night  of  sleep  a  reason¬ 
able  accuracy  of  sleep  stage  classification  is  possible.  However  the  variability 
in  heart  rate  from  night-to-night  for  any  one  individual  produces  unacceptably 
poor  classification  results  on  the  second  night. 
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INTRODUCTION 


Extended  periods  of  sleep  deprivation  commonly  produces  a  decrease  of 
performance  capabilities  at  skilled  tasks  (4,10,13,18)  in  addition  to  un¬ 
healthy  changes  in  personality  profiles  (22,23).  Indeed,  hallucinations 
have  been  observed  in  laboratory  experiments  (9) .  Kleitman  (1 2)  reports 
that  "among  the  effects  of  prolonged  wakefulness  are  irritability  and  mental 
disorganization,  leading  to  daydreaming  and  automatic  behavior,  occasionally 
bordering  on  temporary  insanity."  Decreased  per'crmance  prevents  a  person 
from  meeting  the  requirements  of  many  military  situations  which  require  maxi¬ 
ma?  alertness  and  performance  by  the  on-duty  personnel.  Berry  (4)  reports 
that  fatigue  due  to  inadequate  rest  interf erred  with  the  ability  of  the  astronauts 
to  perform  tasks  in  the  Apollo  VII  and  VIII  missions.  Usually  the  state  of  alert¬ 
ness  and  performance  for  an  individual  is  associated  with  the  amount  of  rest 
and  sleep  he  has  obtained. 

Unfortunately  the  t  .chnical  difficulties  in  obtaining  sleep  information 
have  impeded  sleep  research  outside  the  laboratory.  Classically  sleep  is 
evaluated  from  electroencephalographic  data  (EEG) .  The  instability  of  the  EEG 
electrodes  over  extended  periods  of  time  and  tire  lade  of  an  automated  process 
fer  evaluation  of  the  EEG  has  discouraged  meaningful  research  of  sleep  in 
military  situations. 

One  possible  solution  to  this  problem  would  be  die  development  of  an  alternate 
source  of  sleep  information.  The  source  we  consider  in  this  report  is  beal-te- 
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beat  heart  rate.  This  electrophysiological  measurement  is  more  stable  than 
the  EEG  over  long  periods  of  time  and  there  is  evidence  that  average  heart 
rate  is  influenced  by  sleep.  It  is  the  intent  of  this  report  to  determine  the 
feasibility  of  using  instantaneous  heart  rate  in  an  automated  process  as  an 
indicator  of  the  sleep-wakefulness  cycle. 

Statement  of  Problem 

Pilots  in  flight  are  not  always  able  to  report  accurately  their  physical 
condition,  particularly  with  respect  to  drowsiness-wakefulness,  which 
affects  alertness  and  operational  capability.  A  simple  objective  measure  of 
this  condition  under  operational  conditions  is  desirable.  Any  solution  must 
keep  the  sensory  system  simple  and  must  keep  the  required  transmission  band¬ 
width  small.  One  possibility  is  to  derive  this  wakefulness  information  from 
beat-to-fceat  heart  rate  data,  transmit  the  heart  rate  data  to  ground  stations, 
and  use  computer  analysis  to  determine  sleep  stage  from  heart  rate  derived 
measures . 

In  a  previous  report  from  The  University  of  Texas  at  Austin  (21)  we 
described  the  computation  of  several  different  measures  of  beal-to-beat  heart 
rate.  We  also  computed  the  possible  utility  of  each  of  these  measurements 
to  a  sleep  stage  classification  program.  This  report  examines  computer  classi¬ 
fication  of  sleep  stages  utilizing  the  measures  described  in  our  previous  report. 
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Review  of  Literature 

Since  the  classical  work  of  Dement  and  Kleitman  (8),  the  depth  and  dura¬ 
tion  of  sleep  has  been  determined  by  examination  of  electroencephalographic 
data.  However,  the  recording  of  a  full  night  of  sleep  EEG  on  a  strip  chart 
recorder  results  in  a  bulky  set  of  data.  Further  the  visual  interpretation  of 
these  data  is  a  time-consuming  task  since  an  individual  must  manually  scan 
up  to  1,000  feet  of  strip  chart  record!  Trained  personnel,  using  well-docu¬ 
mented  criteria  to  soore  the  sleep  records,  have  not  produced  a  consistent 
procedure  for  scoring  a  night  of  sleep  with  a  guaranteed  accuracy  better 
than  90%  (24).  Monroe  (16)  reports  inter-rater  consistency  in  scoring  sleep 
records  by  different  specialists  to  be  only  65  percent.  Many  inaccuracies 
may  be  due  to  the  marked  degree  of  subjectivity  that  must  be  used  in 
visually  scanning  lengthy  records.  In  seme  situations  such  as  space  flight 
oandwiath,  weight  considerations  and  poor  electrode  techniques  suggest 
that  an  alternate  signal  to  the  EEG  is  needed. 

Coupling  between  sleep  activity  in  the  brain  and  autonomic  nervous 
system  motor  activity  has  been  documented  by  many  investigators  (5,6,11,15,19). 

Snyder  (20)  reported  significant  changes  (P  <  .05)  in  average  levels  of 
blood  pressures,  respiration  and  heart  rate  between  3tages  1-REM  and 
stage  2.  Data  were  recorded  from  twelve  subjects  for  a  total  of  thirty  nights. 
There  was  a  6%  average  increase  in  heart  rate,  a  7%  average  increase  in  res¬ 
piratory  rate  and  a  4%  average  increase  in  systolic  blood  pressure  from  stage  2 
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to  stage  1-REM.  Significant  changes  (P<  .rj)  d-d  not  occur  between  stage  2 
and  combined  stages  3  and  4. 

Since  die  depth  of  sleep  typically  produces  measurable  changes  in 
autonomic  bodily  functions .  we  anticipated  that  autonomic  activity  could 
be  used  to  describe  depth  of  sleep.  Of  the  organs  under  die  control  of  the 
fiutonomic  nervous  system,  the  heart  has  one  of  the  richest  supplies  of  both 
adrenergic  and  cholinercic  nerve  endings.  The  heart  also  produces  an  elec¬ 
trical  signal  that  is  easy  to  measure  (the  electrocardiogram) .  Most  of  the 
information  supplied  to  the  heart  by  the  autonomic  nervous  system  is  reflected 
in  die  instantaneous  R  to  R  interval  (the  beat-by-beat  heart  rate) .  Only  a 
small  portion  of  the  autonomic  information  supplied  to  the  heart  is  reflected 
in  die  electrocardiographic  wave  shape  Thus  the  beot-by-beat  heart  rate 
should  be  a  useful  measure  in  determining  sleep  stages . 

In  a  study  by  Brooks  (  6  )  six  Individuals  (three  husbands  and  wives) 
were  observed  for  fifty  nights  of  sleep.  Brooks  found  a  10%  average  increase 
in  heart  rate  when  the  depth  of  sleep  lessened  by  one  stage  from  stages  4  to  3 
ot  3  to  2.  He  also  found  a  13.7%  average  increase  in  heart  rate  with  two  stage 
lessening  of  sleep  (4  to  2  and  3  to  1).  A  21 .5%  average  increase  in  heart  rate 
occurred  when  sleep  level  lightened  by  three  stages  (from  4  to  1  or  from  stage  3 
to  wakefulness) .  Brooks  concluded  that  sleep  depth  was  probably  reflected 
mere  in  changes  in  cardiac  cycle  length  U.e.  instantaneous  R  to  R  interval  or 
beat-by-beat  heart  rate)  than  in  the  average  heart  rate  values  he  used. 
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Review  of  Our  Past  Work 

In  our  previous  work,  analysis  of  variance  was  applied  to  several  types 
of  measurements  which  were  obtained  from  instantaneous  heart  rate  data. 
Significant  measures  of  sample  mean  value,  sample  standard  deviation  and 
histogram  were  found  for  both  the  instantaneous  heart  rate  and  the  beat-by¬ 
beat  interval .  Table  1  graphically  illustrates  the  level  of  significance  of 
each  of  these  measures  for  each  subject  fcr  the  interval  measures.  The 
level  cf  significance  is  measured  with  the  F  test  and  indicates  the  level  of 
rejection  of  the  hypotheses  of  equal  mean  value  for  each  stage  of  sleep  for 
the  measure  under  consideration. 

Fourier  analysis  of  the  instantaneous  heart  rate  data  produced  a  large 
number  of  variables  which  were  significant  at  the  .001  level.  All  measures 
were  made  or.  ensembles  of  128  heart  beat  intervals  (the  number  of  seconds 
between  each  heart  beat) .  A  detailed  presentation  of  the  results  of  this  work 
is  available  (21).  At  the  time  we  wrote  this  earlier  report,  it  was  recognized 
that  the  analysis  of  variance  test  did  not  define  the  level  of  separation  a 
variable  might  accomplish  in  multiple  class  data.  However,  variables  that 
are  significant  at  the  .001  level  can  frequently  provide  a  reasonable  starting 
point  in  the  search  for  reliable  measures  to  be  used  for  classification. 


Rationale  of  Present  Approach 


The  data  of  all  test  subjects  presented  to  The  University  of  Texas  was 


classified  into  one  of  six  states  (or  levels)  of  consciousness.  These  were 
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awake-,  stage  one  through  four  sleep,  and  stage  REM.  While  these  categories 
or  stages  of  consciousness  have  been  considered  by  many  to  be  standard,  they 
are  based  primarily  on  EEG  criteria.  Physiologically  these  stages  probably 
have  little  direct  meaning . 

Because  of  the  difficulty  of  computer  classifying  six  categories,  we 
attempted  to  simplify  the  problem.  The  deep  sleep  stages  3  and  4  were  com¬ 
bined  into  a  single  stage  and  stage  l  and  REM  wete  grouped  together.  There¬ 
fore  the  number  of  consciousness  levels  was  reduced  from  6  to  4. 

Stage  1  and  REM  were  combined  because  we  felt  there  was  attle  meaningful 
difference  between  the  two  sleep  stages .  Sleep  stage  REM  is  usually  scored 
whenever  Stage  1  EEG  is  found  after  the  1st  sleep  cycle.  In  other  words,  the 
only  time  Stage  1  was  scored  was  during  the  first  sleep  cycle;  after  that  all 
Stage  1  EEG  was  typically  scored  as  REM.  Agnew  (2)  concurred  in  this  opinion. 

Similarly,  we  were  unable  to  find  any  reason  why  Stages  3  and  4  should  net 
be  combined.  Typically  in  the  evaluation  of  sleep  effectiveness,  stages  3  ;md 
4  are  c  ombined  (2) . 

Once  the  dimensionality  of  our  problem  was  reduced  to  four,  we  sought 
to  develop  a  "cost"  function  to  evaluate  classification  procedure.  This  was 
necessary  since  it  is  typically  difficult  for  humans  to  compare  error  functions 
of  four  variables.  The  weightings  suggested  by  Agnew  (2)  stress  the  importance 
of  differentiating  between  awake  and  sleep. 
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A  much  smaller  cost  is  assigned  for  making  "one"  stage  errors  during  sleep. 
Table  2  represents  the  empirically -derived  cost  matrix. 


TABLE  2 


The  cost  of  classifying  stage - 


Awake 

1,  REM 

2 

3,4  | 

Awake 

0 

1 

1 

i 

1 

i 

t  1  /  hiim 

1 

0 

- j 

1/4  ;  1/2  ' 

_ _  ' _ i 

■ 

2 

1 

1/4 

0 

1/4 

3,4 

1 

1/2 

1/4 

0 

_ 

_ 

A  cost  per  classification  is  obtained  when  this  cost  matrix  is  multiplied  by 
a  sleep  result  matrix  (or  an  element-by-element  basis  -  not  as  matrix  multi¬ 
plication  is  typically  performed)  and  the  elements  of  this  product  matrix  are 
summed  then  divided  by  the  total  number  of  points  considered  in  the  classi¬ 
fication  matrix. 

Note  that  the  greatest  cost  of  errors  occurs  when  an  awake  epoch  is 
classified  as  sleep,  or  when  any  sleep  epoch  is  classified  as  awake.  No 
cost  ii  accrued  for  a  correct  answer  and  weightings  of  one-half  are  charged 
for  missing  the  sleep  stage  by  more  than  one  levei.  This  single  number  cost 
value  permits  effective  comparison  of  the  various  variables  used  in  the  classi¬ 


fication  procedure. 
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DISCRIMINANT  ANALYSIS 


The  optimum  procedure  for  classifying  a  sample  of  data  with  d-measures 
is  achieved  using  the  Bayes  discriminant  function  which  requires  knowledge 
(or  estimates)  of  the  .3  priori  probability  of  occurrence  of  each  category  and 
the  d-dimensional  joint  density  function  of  the  measures  for  each  category. 

A  non -optimal ,  but  computationally  feasible,  approach  to  classification  is 
Multiple  Discriminant  Analysis .  This  analysis  includes  a  linear  transformation 
to  reduce  the  dimensionality  of  the  problem  and  a  Bayes  classifier  as  illus¬ 
trated  in  Figure  1 


Data  X 

Discriminant 

Y 

Bayes  Classifier 

-  -  -  -  _ ^ 

analysis  class 

d-dimensional 

separating 

- >» 

measure  vector 

transform 

Classificatior 

- 


CLASSIFICATION  PROCEDURE 
Figure  l 

Often  information  from  an  experiment  can  be  divided  into  a  sequence  of  con¬ 
secutive  epochs .  For  each  epoch  a  set  of  descriptors  or  measures  is  computed 
that  may  contain  sufficient  information  for  the  classification  of  the  epoch. 

It  is  convenient  to  picture  the  set  of  data  for  each  epoch  as  a  point  in  an 
d-dimensional  space  where  one  point  describes  each  epoch .  Further  assume 
the  classification  associated  with  each  of  these  epochs  is  known;  where 


i  is  equal  to  1 ,  2,  3, . .  .R  and  R  is  the  total  number  of  classes. 

Another  important  representation  of  the  d  descriptors  is  that  of  a 
d-dimensional  vector  \.  Each  epoch  is  represented  by  a  different  vector. 

Thus  we  have  the  picture  of  points  in  a  d -dimensional  space  and  their  corres¬ 
ponding  vector  representation  X. 

The  linear  transformation  of  Figure  1  maximizes  the  distance  between 
centroids  of  each  category  in  the  Y  space  while  holding  the  oveiall  variance 
constant.  The  transformation  reduces  the  dimensionality  to  the  minimum 
number  required  to  compartmentize  the  space  for  the  categories  under  con¬ 
sideration.  That  is,  the  dimensionality  of  Y  is  the  minimum  of  either  (a)  dimen 
sions  of  X,  or  (b)  the  number  of  categories  minus  one. 

The  discriminant  analysis  procedure  assumes  the  joint  density  function 
of  Y  for  each  category  i?  normally  distributed.  Thus,  conditional  probability 
of  group  occurrence  may  be  computed  according  to: 


where 


P(i/Y) 


T7 

i 


f. 

i 


R 

I 

1=1 


i—  1,2, ...R 


is  the  a  priori  probability  of  occurrence  of  class  j 
is  the  joint  density  function  for  class  j  evaluated  at  Y. 


1 


{2rr)d/2  |Z  I* 
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is  the  centroid  of  class  j  in  the  Y  space 
2^  is  the  covariance  matrix  of  Y  for  class  j 

Application  to  Heart  Rate  Data 

The  heart  rate  data  was  divided  into  a  sequent  a  of  !28  beat-co-bsat 
epochs.  For  each  epoch  heart  rate  descriptors  were  computed  that  were  antici 
pated  to  contain  sufficient  information  for  the  determination  of  sleep  stage 
during  that  epoch. 
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PROCEDURE 


Data  Acquisition 

Electroencephalographic  and  electrocardiographic  data  for  this  study 
were  collected  on  FM  magnetic  tape  at  The  University  of  Florida  Sleep  Labo¬ 
ratory  by  W.  H,  Agnew,  Jr.  under  Air  Force  Contract  No.  F41S09-68-C-0G3. 

Ten  subjects  at  The  University  of  Florida  were  selected  on  the  basis  of  good 
physical  and  mental  health  as  determined  by  medical  examination  and  the 
Minnesota  Multiphasic  Personality  Inventory.  Each  subject  spent  at  least 
three  consecutive  nights  in  the  laboratory.  The  first  night  was  used  to  con¬ 
dition  the  subjects  to  the  laboratory  in  order  to  avoid  first  night  effects  (  1  ) . 
Eight  hours  of  sleep  EEG,  EGG,  and  electrocardiograms  were  recorded  on  FM 
magnetic  tape  for  two  nights.  The  method  of  data  collection  and  sleep  scoring 
of  the  EEG  was  reported  in  detail  by  Leonard,  Ford  and  Hughes  of  the  USAF 
School  of  Aerospace  Medicine  (14).  Copies  of  the  FM  magnetic  tapes  containing 
the  sleep  data  were  supplied  to  The  University  of  Texas  Bio- Medical  Engineering 
Laboratory  by  the  U.  S.  Air  Force  School  of  Aerospace  Medicine. 

Data  Reduction 

Double  differentiation  of  the  filtered  electrocardiogram  and  threshold  logic 
units  were  used  to  detect  the  peak  of  each  R-wave.  The  time  in  milliseconds 
between  heart  beats  was  written  on  digital  magnetic  tape .  The  data  were 
grouped  into  records  containing  128  consecutive  beat-to-beat  intervals. 


I'tlfawliiWf ^Sainr  i<  i 
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A  number  of  descriptors  were  then  computed  for  each  record.  Eleven  of 

—  2 

the  descriptors  were  the  mean  value  X,  the  sample  variance  S  ,  and  the  nine 
intervals  histogram  of  standardized  instantaneous  beat-to-beat  intervals  for 
each  record.  Standardized  beat-to-beat  intervals  were  calculated  using 
equation  (1) . 


The  R-R  interval  for  each  beat  was  X^^  and  the  mean  value  and  standard  devi¬ 
ation  of  the  128  R-R  intervals  was  X  and  S,  respectively.  The  standardized 
instantaneous  value  for  the  ith  heart  beat  is  Z^.  The  distribution  of  Z  scores 
calculated  in  this  way  always  has  a  mean  value  of  zero  and  standard  deviation 
of  one.  Other  measurements  for  the  128  R-R  intervals  included  eleven  instan¬ 
taneous  heart  rate  measures  analagous  to  the  interval  measures  and  64  Fourier 
Transform  measures. 


Histogram  intervals  of  one-half  standard  deviation  were  selected  as 
Illustrated  in  Table  3 .  The  number  of  standardized  scores  falling  in  each 
of  the  nine  intervals  were  used  as  the  value  of  the  histogram  descriptors. 

Each  group  of  128  heart  intervals  had  a  different  mean  value,  X  and  sample 
standard  deviation,  S.  The  outer  intervals  of  the  histogram  extended  from 
-  »  to  -  1-3/4  and  +  1-3/4  to  +  » . 

Analysis  of  variance  (7,25)  was  used  to  determine  measure  significance  for 
each  subject.  This  procedure  tested  the  null  hypothesis  for  each  measure. 

In  other  words,  what  is  the  probability  that  a  measure's  mean  value  is  the 


sTjta.  ^ 


TABLE  3 

HEART  INTERVAL  MEASURES 
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Measures 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Description  of  the  Measure 

Sample  Mean  Value 

Sample  Standard  Deviation 
Histogram  Measures  (1/2  o  Intervals) 
z  <  -  13/4  a 

-1  3/4  c  <  z  <  -  1  1/4  a 

-1  1/4  o  <  z  <  -  3/4  o 

-  3/4  cr  <  z  <  -  1/4  o 

-  1/4  a  <  e  <  -  1/4  c 

1/4  c  <  z  <  3/4  a 

2/4  a  <  e  <  1  1/4  a 

l  1/4  a  <  e  <  1  3/4  a 


II 


i  3/4  a  <  z 
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same  for  awake  and  the  five  stages  of  sleep.  Our  prcyram  used  an  F-test  to 
compute  the  probability  that  the  mean  values  of  the  descriptors  for  each  stage 
of  sleep  were  equal.  The  hypothesis  was  rejected  if  P  <  .01 .  The  most  sig¬ 
nificant  measures  were  selected  to  be  used  to  train  the  discriminant  function. 
Detailed  results  of  the  analysis  of  variance  procedure  are  in  our  previous 
report  ( 21). 

Discriminant  Analysis 

After  descriptors  for  each  epoch  of  128  heart  beats  had  been  computed, 
computer  classification  was  performed  in  the  following  manner. 

(1)  Approximately  50%  of  the  first  night's  epochs  were  selected 
as  training  data.  The  exact  number  of  epochs  corresponded  to  the  larger  of 
either  25%  of  that  class's  total  number  of  epochs  for  an  individual's  two  nights 
of  sleep  or  ten  epochs.  Ten  epochs  represented  an  arbitrary  minimum  number 

*f  samples  for  estimation  of  die  d-dimensional  centroid  and  covariance  matrix 
for  each  class.  If  the  first  night  of  sleep  did  not  contain  the  required  number 
of  epochs  then  epochs  were  selected  from  the  second  night  of  sleep.  Each 
epoch  contained  a  complete  set  of  heart  rate,  interval  and  Fourier  transform 
descriptors . 

(2)  Selected  subsets  of  descriptors  from  the  training  set  were  entered 
into  the  discriminant  analysis  program  which  evaluated  the  class  separating 
linear  transformation  X  -*  Y  and  computed  mean  values  and  covariance  matrices 
used  in  the  Baysean  conditional  probability  estimates. 


mss. 


“,aa"‘T>-r‘  iTT-frfiifr*^  — ~  rfrnrf  wit 


16 


(3)  The  accuracy  of  selected  subset  of  descriptors  was  pretested  by 
computing  the  conditional  probability  P(i/X)  for  each  epoch  in  the  training  set 
and  comparing  the  results  to  the  EEG  scored  sleep  stages. 

(4;  The  two  nighcs  of  recorded  sleep  then  were  classified  by  trans¬ 
forming  the  selected  subset  of  descriptors  one  epoch  at  a  time  and  computing 
the  conditional  probability  of  classification  for  each  epoch.  The  computer 
classifications  were  recorded  on  magnetic  tape  for  display  and  accuracy  com¬ 
putations  . 

(5)  The  comouter  classified  sleep  patterns  were  plotted  on  a 
CalComp  plotter  and  the  accuracy  was  determined  by 

(a)  presentation  of  an  error  table  and  computation 
of  the  percent  correctness 

(b)  computation  of  the  correlation  coefficient  between  the 
computer  scored  and  EEG  scored  sleep  records 

(c)  pooling  ihe  data  into  four  classes  (i)  awake,  (ii)  1, 

1-REM,  (iii)  2,  and  (iv)  3,  4  and  evaluating  the  average 
cost  per  epoch  based  upon  cost  function  presented  in 
Table  2,  page  7. 

A  flow  graph  of  the  analysis  procedure  is  shown  in  Figure  2. 


Mb IttBiiattK 
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ANALYSIS  PROCEDURE  FOR  EACH  SUBJECT  -  Figure  2 
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RESULTS 

The  results  of  this  study  are  presented  in  both  tabular  and  graphic 
form.  Tables  4  through  ?.2  present  the  tabular  data.  The  even-numbered 
tables  (4,  6,  8,  10,  etc.)  contain  the  results  of  classifying  the  first  night 
of  recorded  sleep  and  the  odd-numbered  tables  contain  the  results  for  the 
second  night  of  recorded  sleep.  In  most  cases,  the  discriminant  algorithm 
was  trained  on  data  obtained  from  the  first  night  of  recorded  sleep.  Results 
are  presented  for  using  both  balanced  (all  classes  equally  likely)  and  un¬ 
balanced  (actual  frequency  of  occurrence  of  each  category)  a  priori  prob¬ 
abilities.  For  the  specified  sets  of  descriptors  the  tables  contain 

(1)  Accuracy  for  both  six  and  four  category  classifications 

(2)  Average  cost  per  epoch  for  the  night 

(3)  Correlation  coefficient  between  EEG  hand-scored  sleep 
heart  rate  computer- scored  sleep. 

In  these  tables.  Variable  1  is  the  mean  beat-to-beat  interval  of  the 
epoch.  Variable  2  is  the  sample  variance  of  each  epoch,  and  Variable  3 
through  11  are  the  interval  histogram  values.  Table  3,  summarizes 
the  relation  between  variable  number  and  physical  beat-to-beat  interval 
measures.  (For  any  eooch  the  sum  of  the  histogram  measuies  is  128). 

The  variables  noted  in  the  tables  as  "11  Fourier"  represent  11  of  the 
best  Fourier  variables  selected  by  analysis  of  variance  (21).  The  measure 
set  "Histogram  and  Four"  combines  the  8  Interval  Histogram  measures  with  the 
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11  Fourier  measures  * 

The  conditional  probabilities  of  each  sleep  stage  were  calculated 
assuming  each  sleep  stage  to  be  equally  likely  (balanced)  ar.d  using  a  priori 
probabilities  actually  based  on  the  individual  subject.  The  a  priori  prob¬ 
abilities  of  each  of  the  sleep  levels  for  each  subject  is  presented  in  Table  24. 

Five  measures  of  merit  of  the  classification  procedure  are  presented  for 
each  combination  of  measures  used  and  for  each  set  of  a  priori  probabilities. 
The  percent  classification  of  the  first  half  and  the  second  half  of  each  night 
using  six  sleep  categories  are  presented.  In  addition,  the  percent  correct 
classification  using  four  categories  is  presented..  The  weighted  cost  per 
epic  is  listed  as  a  measure  o£  machine  scoring  effectiveness.  The  corre¬ 
lation  coefficient  given  is  a  rough  indication  of  the  sameness  in  shape  between 
the  hand-scored  and  machine-scored  data. 

Tables  5  through  23  and  all  odd-numbered  tables  in  between  represent 
the  results  of  discriminant  classification  using  a  discriminant  function  trained 
on  the  first  night  of  sleep  and  used  to  classify  the  second  night  of  sleep.  The 
same  variables  and  a  priori  probabilities  presented  in  the  preceding  even- 
numbered  table  are  used  in  the  odd-numbered  table.  The  measures  of  merit 
are  tire  same  as  those  used  for  the  preceding  even-numbered  tables.  Occa¬ 
sionally  when  insufficient  samples  were  available  in  the  first  night  of  data,  a 
few  samples  had  to  be  obtained  from  the  second  night  of  sleep  to  train  tha  dis¬ 
criminant  function. 


CLASSIFICATION  OF  THE  FIRST  NIGHT  OF  RECORDED 


*Fllo  1  was  the  first  half  oi  the  first  night  of  data  on  each  subject. 

Most  of  the  training  set  came  from  file  1 . 

'♦File  2  was  the  second  half  of  the  first  night  of  data  on  each  subject. 

Soma  of  the  training  set  came  from  this  data  fite. 

t Correlation  between  entire  night  of  steep  as  machine  scored  against  th3  hand  scoied  night  of  sleop. 


TABLE  5 

CLASSIFICATION  OF  THE  SECOND  NIGHT  OF  RECORDED 
SLEEP  WITH  DISCRIMINANT  ANALYSIS  USING  TRAINING  DATA  FROM  THE  FIRST  NIGHT  OF  SLEEr 
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TABLE  7 

CLASSIFICATION  OF  THE  SECOND  NIGHT  OF  RECORDED 
SLEEP  WITH  DISCRIMINANT  ANALYSIS  USING  TRAINING  DATA  FROM  THE  FIRST  NIGHT  OF  SLEEP 


♦roughly  the  first  half  of  the  second  night  of  data  was  on  file  3  for  each  subject. 

♦the  second  half. of  the  second  night's  sleep  was  on  file  4. 

tcorrelatlon  between  the  entire  night  of  hand  scored  sleep  against  the  entire  night  of  machine  scored  sleep. 


TABLE  8 

CLASSIFICATION  OF  THE  FIRST  NIGHT  OF  RECORDED 
SLEEP  BY  A  DISCRIMINANT  FUNCTION 


••File  2  we •  the  second  hall  of  the  first  night  of  data  on  each  subject. 

Some  of  the  training  set  came  front  this  data  file. 

t Correlation  between  entire  night  of  sleep  as  machine  scored  against  the  hand  scored  night  of  steep. 


TABLE  9 

CLASSIFICATION  OF  THE  SECOND  NIGHT  OF  RECORDED 
SLEEP  WITH  DISCRIMINANT  ANALYSIS  USING  TRAINING  nATA  FROM  THE  FIRST  NIGHT  OF  SLEEP 
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fcorrelalion  between  the  entire  night  of  hand  scored  sleep  agalns 


TABLE  13 

CLASSIFICATION  OF  THE  SECOND  NIGHT  OF  RECORDED 
SLEEP  WITH  .DISCRIMINANT  ANALYSIS  USING  TRAINING  DATA  FROM  THE  FIRST  NIGHT  OF  SLEEP 


Most  of  the  training  set  came  from  file  1 . 

**File  2  was  the  second  half  of  the  first  night  of  data  on  each  subject. 

Some  of  the  training  v.et  came  from  this  data  file. 

■f Correlation  between  entire  night  of  sleep  as  machine  scored  against  the  hand  scored  night  of  sleep. 


TABLE  15 

CLASSIFICATION  OF  THE  SECOND  NIGHT  OT  RECORDED 
SLEEP  WITH  DISCRIMINANT  ANALYSIS  USING  TRAINING  DATA  FROM  THF.  FIRST  NIGHT  OF  SLEEP 


tcorrolatipn  botween  the  entire  night  of  hand  scored  sleep  against  the  entire  night  of  machine  scored  sleep. 


♦Fite  1  was  the  first  half  of  the  first  night  of  data  on  each  subject. 

Most  of  the  training  set  came  from  file  1. 

**FUe  2  was  the  second  halt  of  the  first  night  of  data  on  each  subject. 

Some  of  the  training  set  came  from  this  data  file. 

tCorrelation  between  entire  night  of  steep  as  machine  scored  against  the  hand  scored  night  of  sleep. 


CLASSIFICATION  OF  THE  SECOND  NIGHT  OF  RECOPDED 
SLEEP  WITH  DISCRIMINANT  ANALYSIS  USING  TRAINING  DATA  FROM  THE  FIRST  NIGHT  OF  SLEEP 


CLASSIFICATION  OF  THE  FIRST  NIGHT  OF  RECORDED 
SLEEP  BY  A  DISCRIMINANT  FUNCTION 


•File  l  was  the  first  half  of  the  first  night  of  data  on  each  subject. 

Most  of  the  training  set  came  from  file  1 . 

•File  2  v.-as  the  second  half  of  the  first  night  of  data  on  each  subject. 

Some  of  the  training  set  came  from  this  data  file. 

^Correlation  between  entire  night  of  sleep  as  machine  scored  against  the  hand  scored  night  of  steep. 


CLASSIFICATION  OF  THL  SLCOND  NIGHT  OF  RECORDED 
SLEEP  BY  A  DISCRIMINANT  FUNCTION 


TABLE  22 

LASSIFICATION  OF  THE  FIRST  NIGHT  OF  RECORDED 
SLEEP  Btf  A  DISCRIMINANT  FUNCTION 
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Most  of  the  training  set  came  from  file  I . 

•♦File  2  was  the  second  half  of  the  first  night  of  data  on  each  subject. 

Some  of  the  training  set  came  from  this  data  file. 

tCorrelation  between  entire  night  of  sleep  as  machine  scored  against  the  hand  scored  night  of  sleep. 


TABLE  23 

CLASSIFICATION  OF  THE  SECOND  NIGHT  OF  RECORDED 
SLEEP  WITH  DISCRIMINANT  ANALYSIS  USING  TRAINING  DATA  FROM  THE  FIRST  NIGHT  OF  SLEEP 


tcomalatlon  between  the  entire  night  of  hand  scored  ■  nst  the  entire  night  of  machine  scored  sleep. 
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TABLE  24 

PROBABILITY  OF  OCCURRENCE,  P 


A 

1 

2 

3 

4 

REM 

Safer 

.237 

.072 

.328 

.072 

.072 

.208 

Schmidt 

.090 

.090 

.423 

.090 

.090 

.208 

Farrington 

.087 

.087 

.434 

.087 

.087 

.217 

Chinoy 

.088 

.088 

.377 

.079 

.088 

.281 

Gildersleeve 

.082 

.002 

.418 

.082 

.098 

.237 

Moss 

.095 

.103 

.465 

.086 

.000 

250 

Verick 

.076 

.076 

.481 

.076 

.091 

•  )8 

Nordyke 

.088 

.088 

.368 

.088 

.158 

.210 

Phillips 

.109 

.099 

.446 

.099 

.093 

.148 

Padula 

.075 

.090 

.466 

.075 

.075 

.218 

Number  of  Epochs  of  Stage  I 

P  =  In  Training  Set _ 

Total  Number  of 
Training  Epochs 

I  =  Awake,  1,  2,  3,  4,  and  REM 
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Tables  25,  26,  27  and  28  represent  examples  of  good  results  using  the 
discriminant  function  procedure. 

Tables  25  and  26  represent  results  obtained  from  an  attempt  to  classify 
the  first  and  second  nights  of  subject  Chinoy's  sleep  with  variables  1,2, 3, 5, 7, 
9  and  11.  A  graphical  representation  of  his  sleep  pattern  for  the  two  nights 
based  on  EEG  hand-scored  records  is  presented  in  Figures  3  and  4.  The  corres¬ 
ponding  heart  rate  machine  scored  sleep  patterns  are  shown  in  Figures  5  and  6. 

Tables  27  and  28  give  the  details  of  an  analysis  of  the  same  data  using 
a  combined  histogram  Fourier  analysis  measure  set.  The  computer-generated 
classification  based  on  these  measures  is  presented  in  Figure  7  (first  night  of 
recorded  sleep)  and  Figure  8  (second  night  of  recorded  sleep) . 


— -c~--  n^-  .t.Yi-ii  y-i 
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TABLE  25 


EXAMPLE  OF  CLASSIFICATION  RESULTS  FOk  FIRST  NIGHT  OF  RECORDED  SLEEP 


Subject: 

Variables: 

Night: 

A  priori  probability: 


Chinoy 

1,2,3,5,7,9,11 
3  (first  recorded  night  of  sleep) 
Actual  frequency  of  occurrence 


This  stage  of  sleep. 


was 

classified  as  Awake 


Awake 


REDUCED  SLEEP  MATRIX 


This  stage  of  sleep. 


Awake 


(1  ,R£M) 


(3,4) 


was 

classified 


Awake 
(1  ,REM) 


(3,4) 


j'lfidtt&.U  gg 


13 


TABLE  26 

EXAMPLE  OF  CLASSIFICATION  RESULTS  FOR  SECOND  NIGHT  OF  RECORDED  SLEEP 


S'vtbject: 

Variables: 

Night: 

A  priori  probability: 


Chinoy 

1,2,3,5,7,9,11 

4  (second  recorded  night  of  sleep) 
Actual  frequency  of  occurrence 


• 

This  stage  of  sleep. , 

• 

• 

Awake 

1 

2 

3 

4 

REM 

was 

classified 

Awake 

5 

0 

1 

0 

0 

5 

as 

1 

1 

1 

3 

0 

G 

10 

2 

0 

2 

73 

3 

14 

16 

3 

1 

1 

4 

1 

2 

3 

4 

0 

0 

7 

0 

0 

0 

REM 

0 

1 

1 

0 

0 

27 

REDUCED  SLEEP  MATRIX 

This  stage  of  sleep. 

• 

• 

Awake 

(1 ,  REM) 

2 

(3,4) 

• 

was 

Awake 

5 

5 

1 

0 

classified 

as 

(I  ,REM) 

1 

39 

A 

0 

2 

0 

18 

78 

17 

(3,4) 

1 

4 

11 

3 

ittsM-iigifaft  t  'itrt-  1  ‘i  i~i  i 


Night  '  3 
Tapes  :  2S,  30 
Start  Time:  0849 
End  Time:  1647 


\  !  \  t  nr  nr  iTnrnTr  :  tt i  r~r  ft 


Time 

Tic  marks  =  15  min 

Above  Stage  i 
indicates  S.E.M. 


Marks  points  of 
a  ^agreement  on  stage 


E£6 

SLEEP  PATTERNS 

Figure  3 


CHINOY 
Night  •  4 
Top«s  :  35,32 
Storl  Tint: 0857 
End  Tim* :  1654 


J 


Mark*  points  of 

Tim* 

disogrtomon!  on  stag* 

Tic  marks  =15  min 

Abo**  Stoge  I 
indicates  RJE.M. 

EEG 

SLEEP  PATTERNS 

Figure  4 

Above  Stage  ! 
indicates  RE.M. 


Time  Tic  Marks  *  15  win 


CLASSIFICATION  OF  HEART  RATE  DATA 
INTO  SLEEP  STAGES 


Figure  5 


CHINOY 
NIGHT  4 


VARIABLES  1,2,3.5,7,9,11 


firm  i  i  m  i  i  ii  ri  i  ttt  rrn  r  nr 

- Above  Stage  I  Time  Tic  Morks  »  15  min 

intiieotes  R.  EM. 


CLASSIFICATION  OF  HEART  RATE  DATA 
INTO  SLEEP  STAGES 

Figure  6 
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TABLE  27 

EXAMPLE  OF  CLASSIFICATION  RESULTS  FOR  FIRST  NIGHT  OF  RECORDED  SLEEP 


Subject:  Chinoy 

Variables:  1,  4-11  and  eleven  Fourier 

Night:  3 

A  priori  probability:  Actual  frequency  of  occurrence 

This  stage  of  sleep . 


• 

Awake 

1 

2 

3 

4 

REM 

was 

classified 

Awake 

4 

2 

2 

1 

0 

1 

as 

1 

0 

3 

5 

0 

0 

3 

2 

0 

4 

55 

3 

14 

6 

3 

0 

0 

3 

1 

1 

0 

4 

1 

0 

5 

0 

9 

0 

REM 

2 

4 

5 

0 

0 

58 

REDUCED  SLEEP  MATRIX 

This  stage  of  sleep. 

• 

• 

Awake 

(1,  REM) 

2 

(3,4) 

* 

was 

Awake 

4 

3 

2 

1 

classified 

as 

(1  ,REM) 

2 

68 

10 

0 

2 

0 

10 

55 

17 

(3,4) 

1 

0 

8 

11 
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TABLE  28 


EXAMPLE  OF  CLASSIFICATION  RESULTS  FOR  SECOND  NIGHT  OF  RECORDED  SLEEP 


Subject: 

Variables: 

Night: 

A  priori  probability: 


Chinoy 

1,  4-11  and  eleven  Fourier 

4 

Actual  frequency  of  occurrence 


This  stage  of  sleep 


Awake  12  3 


was  Awake 

classified 

as 

2 

3 


0 

1 

1 

0 


0  1  1 

1  9  0 

4  74  3 

0  2  0 


4  0 

REM  C 


0  6  0 
0  1  0 


4 

0 

1 

12 

2 

1 

0 


REM 

3 

11 

17 

0 

0 

30 


REDUCED  SLEEP  MATRIX 

This  stage  of  sleep . 

Awake  (l,  REM)  2  (3,4) 


was 

Awake 

0 

3 

1 

1 

classified 

as 

(1  ,REM) 

l 

42 

10 

1 

2 

1 

21 

74 

15 

(3,4) 

0 

0 

8 

3 

CH  INOY 
NIGHT  3 


VARIABLES:  HISTOGRAM  AND 

FOURIER 


t-H-H-H  I  H  +4-H  I'HMIIHI  t-H-H-H 


Above  Stogc  l  Time  Tic  Marks  «  15  min 

indicates  R.E.M 


CLASSIFICATION  OF  HEART  RATE  DATA 
INTO  SLEEP  STAGES 

Figure  7 


CHiNOY 
NIGHT  4 


VARIABLES:  HISTOGRAM  AND 

FOURIER 


H-H  K  H  W’t+4-H-H  M  MH  W4-f4-H-4-4 


Above  Stage  l  Time  Tic  Marks  *  15  min. 

indicates  R.E.M. 


CLASSIFICATION  OF  HEART  RATE  DATA 
INTO  SLEEP  STAGES 


Figure  8 


52 


DISCUSSION 

The  overall  results  of  a  classification  by  discriminant  analysis 
were  somewhat  discouraging  when  the  second  night  of  data  was  considered. 
As  might  be  expected  the  unbalanced  a  priori  probabilities  consistenUy 
produced  better  classification  results. 

Much  of  the  inability  of  this  algorithm  to  correctly  classify  sleep 
is  rooted  in  the  fact  that  there  was  considerable  amount  of  intra-subject 
variation  and  even  intra-night  variation  in  the  beat-to-beat  heart  rate. 
Aldredge ,  et  al  (  3)  has  investigated  the  intra-subject  and  intra-cycle  vari¬ 
ation  in  the  mean  heart  rate  and  sample  standard  deviation  of  the  data  used 
for  this  investigation.  They  examined  the  possibility  that  the  mean  values 
of  average  heart  rate  and  sample  standard  deviation  might  not  be  consistent 
throughout  a  night  of  sleep  or  between  two  nights  of  sleep.  Undesired 
variations  in  the  average  and  sample  standard  deviation  in  heart  rate  during 
a  night  of  sleep  were  determined  by  testing  the  following  hypothesis  with 
analysis  of  variances  for  each  subject: 

H  :  The  mean  values  of  a  random  variable  X  for  each 
o 

stage  I  are  equal  for  cycles  of  sleep  during 
a  single  night  of  sleep. 

{ X  is  either  equal  average  heart  rate  or  sample  standard 
deviation,  and  I  represents  either  (X,REM),  2,  or  {3, 4))» 
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Aldredge,  et  al,  concluded  that  for  a  single  night  of  sleep  the 
above  hypothesis  could  be  rejected  in  most  of  the  cases  at  .0001  sig¬ 
nificance  level  when  X  represented  the  average  heart  rate  and  I  was  either 
equal  to  (i,REM)  or  2  or  (3,4).  However,  the  hypothesis  could  not  be 
rejected  at  the  .0001  level  for  any  stage  of  sleep  when  the  sample  stan¬ 
dard  deviation  was  tested.  A  close  inspection  of  the  intr-cycle  variation 
in  average  heart  rate  suggested  that  the  heart  rate  during  a  cycle  of  sleep 
was  influenced  by  the  average  REM  heart  rate  at  the  onset  of  the  cycle. 

The  null  hypothesis  of  equal  means  for  average  heart  rate  values  of  stage 
REM  and  of  stage  2  for  any  given  cycle  was  rejected  in  favor  of  the  alter¬ 
nate  hypothesis  that  the  mean  for  REM  was  greater  than  the  mean  for  stage 
2  at  the  .05  significance  level.  Also  the  alternate  hypothesis  was  accepted 
when  the  mean  values  of  mean  averaged  heart  rate  values  of  stage  REM  are 
greater  than  that  of  combined  stage  REM  and  combined  stages  3  and  4  were 
compared  with  a  one  tailed  t-test  at  a  5%  significance  level.  A  third  null 
hypothesis  which  stated  that  the  mean  averaged  heart  rate  values  of  stage  2 
were  equal  to  that  of  combined  stages  3  and  4  could  not  be  rejected. 

The  performance  of  the  Fourier  measure  and  their  ability  to  represent 
sleep  information  was  disappointing.  With  the  except  -tj  of  Messrs  ,  Paduia 
and  Moss,  the  Fourier  measures  alone  did  not  improve  the  cost  of  sleep 
scoring.  When  combined  with  the  histogram  measures  ,  five  of  the  sub¬ 
jects  had  a  slight  decrease  in  cost  of  classification.  Table  29  illu¬ 
strates  the  relative  cost  of  of  using  histogram  measures,  Fourier  measures 


Subject 

Histogram 

Measures 

Fourier 

Measures 

Combined  Fourier  and 
Histogram  Measures 

Chinoy 

.12 

.25 

11 

Safer 

.39 

.50 

.40 

Gildersleeve 

.15 

.20 

.17 

Padula 

.50 

.27 

.48 

Moss 

.35 

.19 

.30 

Phillips 

.21 

.23 

.20 

Schmidt 

.25 

.31 

.27 

Farrington 

.14 

.25 

.13 

jS 

I 

1 

5| 

’? 


3 


/. 

i 

4 

4 

% 

•3 

\ 

1 

i 

i 
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alone  and  the  combined  histogram  and  Fourier  measures. 

Examination  of  the  summary  tables  indicated  that  all  nine  histograms 
were  never  used  in  the  same  covaraince  matrix  since  the  linear  dependence 
of  one  histogram  variable  on  the  other  histogram  variables  produces  a  sin¬ 
gular  matrix  which  is  non-invertible. 

Table  30  considers  that  advantages  of  using  balanced  and  unbalanced 
a, priori  probabilities.  In  terms  of  our  empirically-derived  cost  function,  the 
non-balanced  a  priori  probabilities  are  quite  similar.  If  percent  of  accuracy 
diagnosed  sleep  epochs  is  the  criterion,  the  non-balanced  a  priori  probabilities 
are  decidedly  superior. 

CONCLUSIONS 

Result  of  this  research  indicates  that  it  is  possible  to  classify  heart 
rate  patterns  into  sleep  stages.  However,  the  results  are  not  overwhelming, 
in  spite  of  the  fact  that  the  analysis  of  variance  indicates  an  optimistic  possi¬ 
bility  of  sleep  stage  classification  ability  of  these  measures.  Much  of  the 
difficulty  experienced  by  this  algorithm  can  be  attributed  to  Intra-night,  intra¬ 
subject  variations  in  mean  value  as  the  90-minutes  sleep  cycles  progress  through^ 
out  the  night. 

Our  empirically  derived  cost  function  proved  to  be  a  useful  measure  of 
the  effectiveness  of  our  sleep  classification  algorithm.  We  contend  that  any 
measure  of  merit  used  to  evaluate  a  device  or  procedure  should  of  necessity  be 
closely  related  to  the  original  problem,  in  this  case  the  study  of  sleeping  patterns. 
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TABLE  30 

EFFECT  OF  USING  A  PRIORI  PROBABILITIES 
IN  CLASSIFICATION  ALGORITHM 

Classification  Criteria  Number  of  ClassHi cation  I  ans 

I  Lowest  Cost  Using  Night  1  Night  2 

A  Priori  Probability  14  14 

Balanced  Probabilities  0  2 

(likelihood  ratio) 

No  Significant  Difference  4  4 

Between  Above  Methods 

II  Highest  Correlation  Using 


A  Priori  Probability  8  4 

Likelihood  Ratio  7  7 

No  Difference  4  8 

III  Percent  Accuracy  Using 

A  Priori  Probability  15  13 

Likelihood  Ratio  0  0 


No  Difference 


3 


5 


We  also  conclude  that  if  it  were  possible  to  normalized  the  heart  rate 
data  so  that  the  inter  (90  minutes)  sleep  cycle  vailation  in  mean  were  less, 
then  better  results  would  have  been  obtained.  We  suggest  that  if  ar.  alternate 
algorithm  could  be  developed  to  determine  the  beginning  of  each  90  minute 
sleep  cycle,  then  this  algorithm  would  be  able  to  accurately  classify  the 
data  into  sleep  stages. 

Also,  if  it  were  possible  to  extract  respiration  information  from  ampli¬ 
tude  variations  in  the  QRS  complex,  this  additional  variable  might  improve  the 


accuracy  of  the  algorithm. 


Si  igjgSgfNiMiC-j, 
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Modif'cd  Discriminant  Analysis 

The  original  discriminant  analysis  program  is  that  of  Dr.  Donald 
J.  Veldman  of  The  University  of  Texas  at  Austin  ,  Department  of  Educational 
Psychology  and  is  well  documented  in  his  book  Fortran  Programming  for  the 
Behavioral  Sciences.  In  addition  to  his  discriminant  analysis  program,  we 
have  found  it  helpful  to  use  some  of  the  minor  subroutines  included  in  his 
work;  namely  PRTS  and  PCDS,  which  are  efficient  print  and  punch  sub¬ 
routines  and  COR S  which  computes  means,  sigmas,  and  intercorrelations. 
The  discriminant  program  determines  a  transformation  matrix  of  data  from 
known  groups  (sleep  stages)  which  maximizes  the  distance  between  centroid 
while  holding  the  overall  distance  between  points  of  data  constant.  The 
dimensionality  of  the  new  space  is  the  minimum  of  (1)  the  number  of 
groups  minus  one,  or  (2)  the  number  of  variables.  Basically  discriminant 
analysis  is  a  transformation  from  the  data  space  to  a  new  reduced  measure 
space . 

rhe  classification  program  of  Cooley  and  Lohnes  requires  not  only 
the  D-WTS  matrix  and  centroids,  but  the  covariance  matrix  at  each  group  . 
The  modifications  of  the  discriminant  analysis  program  for  computing  the 
covariance  matrices  are  described  in  the  following  paragraphs. 

In  loop  35  of  the  program  original  SDSCR1M  variable  sums  and 
cross  products  are  accumulated  fcr  all  samples  of  the  variables  (i.e.  all 
data  pooled  regardless  of  group  classification) .  After  storing  the  sum  of 
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the  score  for  each  group  in  S  and  the  sum  of  the  cross  products  in  A,  a 
"within  group  matrix"  is  computed  and  added  to  W.  It  is  this  "within 
group  matrix"  for  each  group  which  we  desire  to  use  to  calculate  the  dis¬ 
persion  matrix  (i.e.  covariance  matrix)  for  each  group.  We  wish  to  save 
that  matrix  for  each  group  before  it  is  added  with  those  of  the  other  groups 
into  matrix  W.  This  is  done  in  the  modified  version  called  JDSCRIM  in 
line  c  by  saving  matrix  B  in  CC  which  will  be  added  in  to  W  in  statement 
34.  The  triangular  matrix  we  wish  is  now  available  in  B(I.J.);  we  need  to 
divide  it  again  by  C(M)  [statement  d]  and  fill  in  the  portion  below  the 
diagonal  [statement  e].  If  we  wish  to  calculate  and  punch  the  dispersion 
matrix  for  each  group  a  minus  one  is  placed  in  cc  24-25  of  the  control  card. 
The  matrix  B  for  each  group  is  written  onto  scratch  tape  1  in  loop  33.  After 
the  D-WTS  are  calculated  (through  statement  60)  we  have  the  items  necessary 
for  calculating  the  reduced  space  dispersion  matrices  which  will  be  done  in 
the  subroutine  DMG,  called  after  statement  60. 
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PROGR AM  JO ! SC ( 1 NPUT # OUTPUT/ PUNCH/ TAPE l / TAPES ) 

CALL  JOSCRIM 
CALL  EXIT  *  ENO 
SUBROUTINE  JDSCRIM 

ODIMEnSION  A(7Q#70)#  W(70*70W  C<70/70);  S(70#25?*  T(7Q>/ 

1  V(70>/  X{ 70)/  Y(70)z  Z(70)/  Q(70)/  G(25).  KF(X6)/  KH(X5> 
OIHENSION  8 ( 70/ 70 ) / CC ( 70/ 70 ) 

NX  *  70  *  NS  ■  25 

5  CALL  CCDS  (KF#  NV/  NO,  KW,  KJ,  KEY) 

iP(key,eo..h  rewind  j  ' 

CALL  INPUT  CD/  X,  0/  K F#  NV) 

00  10  I  •  X/NV 
00  XO  J  •  1/NV 

ec/j)  «  o.o 

XO  W(I,J)  «  0.0 

IF  <<T  .OT.  0)  REWIND  2 
00  35  M  «  i/NQ 
READ  15/  N/  KH 
X5  FORMAT  (15/  15A5) 

PRINT  20/  H/  N/  KH 

20  FORMAT  (/  6H  GROUP#  12/  18/  10H  SUBJECTS*/  2X/  X5A5) 

0(M)  *  N 
00  25  I  •  X/NV 
SC/M)  «  0.0 
DO  25  J  •  J/NV 
25  AC/J)  •  0.0 
DO  30  I  •  1*N 

call  input  (ID/  */  n  ♦  m  •  xooo#  KFz  nv) 

IF  (KT  .GT.  0)  WRITE  (£)  10/  (X(j)z  J  •  l/NV) 

DO  30  J  •  1«NV 
S(j/M)  •  S { J/M)  ♦  X(J) 

DO  30  K  *  JjNV 

30  A(J/K)  •  A ( J# K )  ♦  X(J)  •  X ( K ) 

DO  34  I»i/NV 
DO  34  J»IzNV 
CC/J)  ■  CC/J)  ♦  AC/J) 

BC/J)*  (AC/J)  •  SC/M)  •  S ( J/M)  /G(M) ) 

CC(t,J)»  B( 1/ J) 

BC/J)  •  8( !/ J)/3(M! 

B( J#  X )  «B(I/J) 

34  W(I,j)  •  W(I#J)  ♦  CCC/J) 

IF  (KEY.NE.«;»  GO  TO  35 
DO  33  I*X#NV 

33  WRITE  C)  (B(1#J)/J«X/NV) 

35  CONTINUE 

TN  v  SUHF(G/  X#  NO#  N2) 

DO  40  I  •  i#NV 

Tin  *  SUMFIS#  •!/  NG/  NX)  /  TN 
40  0(1)  •  CC/I) 

06  45  I  •  1*NV 
DO  45  J  •  I*NV 

CC/J)  *  CC/J)  /  7N  •  Till  •  Tl J) 

CIJ#D  *  CC#J) 

AC,J)  «  CC/J)  •  TN  •  WC/J) 

4(J/D  «  A(I/J) 

45  W(J#I)  «  W(I/J) 

CALL  INVS  (NV#  W#  X,  Y,  l,  N1 ) 

DO  55  I  •  X#NV 
DO  50  J  •  1#NV 
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: 

l 

C 


f 

I 


50  xtj)  •  w(i# j) 

08  55  J  •  1/NV 

55  W(J,J)  •  SCPF(X,  A/  1#  Uz  NV,  Nl) 

NF  •  MJN0(NG  •  1/  MV) 

CALL  AEVS  (NV,  NF,  0*0/  W,  A /  V/  X/  V/  Z,  Nl) 

03  60  J  •  l/NF 
E  *  1*0  /  S0RT(V(J)) 

09  0  I  •  i/NV 

60  A(!/J)  '  A ( I, J)  •  E 

IF  { <H  .EQ.  1)  CALL  PCDS  (Az  NVz  NF*  5H0  WTSz  NJ) 

IF  (KH  «E0*  1)  CALL  PRTS  (Az  NVz  NFz  5M0  WTSz  Nl) 

IF(KeV«EQ*«1)CALL  DMG(A,NV,nF,nG,NI) 

09  65  I  •  1/iW 
65  X(I)  •  SQi?T(C(IzI)  > 

CALL  AXBS  (Cz  A,  W,  NVz  NF,  NV,  Nl) 

09  70  I  «  1/NF 

70  Y< I )  «  S0RT(SCPF(Az  «z  Iz  I z  NVz  Nl)) 

D8  75  I  •  1/NV 
08  75  J  ■  l/NF 

75  C(),J)  •  W(IzJ)  /  (X(I)  •  Y(J>) 

TR  *  SuMF(Vz  tz  NF,  Nl) 

XL  -  1.0 
09  *0  I  •  lzNF 

xii)  ■  v«n  /  tr  •  ioo.o 

so  xl  «  xl  •  d*o  /  (i*o  ♦  vm>)  »vn*nv»gn-nq*gm«gn«i*o 

SS  •  SORT! ( VN»«2  *  QM*«2  •  4*0)  /  (VN»*2  ♦  GM«m»2  •  5*0>) 

YY  *  XL**(1*0  /  SS) 

FA  *  Yn  •  QM 

FB  »  ((TN  -  1*0)  •  (VN  ♦  ON)  /  2*0)  •  SS  •  (VN  •  GM  •  2.0)  /  2*0 
F  •  <FB  *  (1*0  •  YY))  /  (YY  •  PA) 

P  -  PRBF(FAz  FB*  F> 

PRINT  S5z  XLz  F Az  FBz  Fz  P 

850F8RMaT  (//  I5H  WILKS  LAMBDA  ■/  F10.3  //  7H  D.F,  •  ,  F5,0z 
14H  AND,  F7.0  //  10H  F.RATIB  •*  FS.3z  5X*  3«P  •  *  F7,*) 

OF  *  VN  ♦  ON 

CC  ■  TN  •  OF  /  2.0 

08  90  I  •  l*NF 

CS  •  CC  •  aL8G( 1*0  ♦  V( I ) ) 

DF  •  OF  •  2*0 

P  -  PRBF (OF*  1000. Oz  CS  /  OF) 

90  PRINT  95,  Iz  X(l),  CS/  DF,  P 

950F9RMAT  (/  5HOR88T,  12,  FiO.2,  14H  PCI.  VARIANCE  // 

113H  cHI'SQUARE  f,  F10.3/  5X,  6H0.F.  •«  FS.Oz  5X,  3HP  *z  F?.*> 

08  100  I  •  tzNV 
TU)  «  TO  •  TN 
09  106  <J  •  1/KQ 
100  S(I.J)  •  S(IzJ)  /  G( J) 

CALL  AXBS  (S,  A,  W,  *NCi  NF,  NV,  Nl) 

CALL  PRTS  '.W,  NGz  NF,  5HCENT.,  NJ) 

CALL  PCOS  (Wz  NG*  NF,  5HCENT.,  Nl) 

CALL  PRTS  (Cz  NV,  NFz  6HC8REL.,  Nl) 

OFW  1  TN  *  ON 
PRINT  105,  GM,  OFW 

1050F8RMaT  (//  26H  UNIVARIATE  F-TESTS.  OFB  •,  F3.0z 
16H  OF*  »*  F6.C  ///  IBM  VARIABLE  F.RATJ8,  6X,  IMP) 

08  115  I  •  IzNV 
B  •  0*0 

09  HO  J  •  1/NQ 
110  8  •  B  ♦  S( I* J)**2  *  G(J) 

CC  •  T(I)*»2  /  TN 


■AMHkaHiiiaflk 
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SI 


F  •  <(B  -  CC)  •  DFW)  **  UOtl)  •  B)  •  GH) 
P  •  PRBF<GM#  DFW#  F) 

U5  PRINT  120#  I#  F#  P 

120  FORMAT  {/  16#  F12.A,  FJO#*) 

CALU  PRTS  (S#  NV#  NG#  6HG  MEAN#  NX  > 

IF  (KT  ,£Q,  0)  GO  TO  5 

REWIND  2  *  NT  •  TN 

DO  130  I  •  l#NT 

READ  (2)  ID#  (X{J)#  J  •  1#NV) 

DO  125  J  •  1#NF 

125  Y(J)  •  SCPFJX#  A#  ip  J#  NV,  Nl) 

130  CALI  SyBS  (Y#  NF#  2HDS#  ID) 

GO  TO  5  §  RETURN 

END 
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SUBROUTINE  SDJSCRJM 

00IHENSI8N  A(70,70),  W(70,70),  C(70,70),  S(7Q,25),  T(70i# 
l  V(70),  X(70),  Y<70>,  Z(70),  Q(70),  6(25),  KF (16),  XH(15) 
DIMENSION  B(7C,70),CC(?Q,70) 

N1  p  70  *  N2  ■  25 

5  CALL  CCDS  (KF,  NV,  N0,  KW,  KT,  KEY) 

CALL  INPUT  (ID,  X,  o*  KF,  NV) 

08  10  I  «  l,NV 
08  10  J  •  J,NV 
C(I,J)  •  0.0 
10  WU,j)  •  0.0 

IF  (KT  ,QT,  0)  REWIND  2 
08  35  H  •  l,NQ 
READ  15,  N,  KH 
15  FORMAT  (15,  15A5) 

PRINT  20,  M,  N,  KH 

20  FORMAT  (/  6H  GROUP,  !2,  18,  10H  SUBJECTS.,  2X,  j5A5) 

Q(M)  •  N 
08  25  I  •  i,NV 
S( I,M)  •  0.0 
D8  25  J  •  I,NV 
25  A( I,  j)  •  0.0 
08  30  I  •  1,N 

call  input  do,  x,  n  ♦  m  •  1000,  kf,  nv) 

IF  (KT  .OT.  0)  WRITE  (2)  10,  (X(J),  J  •  l,NV> 

08  30  J  •  1,NV 
S(d,M)  »  S(J,M)  ♦  X(d) 

08  30  K  •  J, NV 

3C  A(J,K)  •  A { J,K )  ♦  X(J)  •  X(K) 

08  3?  I  •  l,NV 
08  35  d  •  I,NV 
Cd,j)  *  C(I,J>  ♦  A ( I , d ) 

35  W(I,J)  ^  W(I,J>  ♦  (A(I,d)  .  S(I,M)  •  S(d,H)  /  Q(M) ) 

TN  •  SUMF(G,  1,  NG,  N2) 

08  AO  I  •  l,NV 

T( I )  •  SUHF(S,  -I,  NO,  Nl)  /  TN 
AO  Q( I )  i  C(1,I) 

08  A5  I  ■  1,NV 
08  A5  d  •  I,NV 

C(I,J)  •  CU,d)  /  TN  •  T(  I )  *  T (d) 

C(d, I >  •  C(I,J> 

A(I,J)  *  C(I,d>  •  TN  .  W(!,J) 

A(J,  J  )  •  Ad,J) 

A5  W(d,I)  •  W(I,d) 

CALL  INV5  (NV,  W,  X,  Y,  2,  Ni ) 

08  55  I  •  1#NV 
08  50  d  •  t,NV 
50  X(d)  «  W(!,J) 

08  55  d  ■  l,NV 

55  W(I,dI  •  SCPFlX,  A,  1,  d,  NV,  Nl) 

NF  o  MinQ(nG  •  1,  NV) 

CALL  AEVS  (NV,  NF,  0#0,  W,  A,  V,  X,  Y,  2,  Nl) 

08  60  d  •  l*NF 
E  •  1*0  /  SQRT(V(d)) 

08  60  I  •  1«NV 
60  AJI,J>  •  A( I,  J)  •  E 

IF  (KW  ,E0.  1)  CALL  PCOS  (A,  NV,  NF,  5HD  WTS,  Nl) 

IF  (KW  «EQ.  1)  CALL  PRTS  (A,  NV,  NF,  5H0  WTS,  Nl) 

08  65*1  •  l#NV 


65  X(I)  •  SQRT<C(t#I) ) 

CALL,  AxBS  (C,  A#  W#  NV#  NF#  NV,  Ni) 

08  70  I  •  t#NP 

70  YU)  •  SQRUSCPFtA#  W#  I#  I#  NV/  Ni)) 
08  75  l  •  1#NV 
08  75  J  •  1#NF 

75  CU#J)  •  W(I/J)  /  (XU)  •  Y(J)) 

TR  i  SU*F(V#  1/  NF#  Nl ) 

XL  *  i.O 
08  80  I  *  1#NF 
XU)  •  VC  I )  /  TR 

#0  XL  •  XL  *  U»0  / 

SQRf  ( (VN«»2 


SS 

YY 

Fa 

FB 

F 

P 


•  100*0 
(1*0  ♦  VU)>) 

•  QM**2  •  A*0> 


♦VN9NV6QN-NQ8QM-3N-1 .0 
/  (VN®*2  ♦  6H**2’»  5*0)) 


XL**(1*0  /  SS) 

VN  •  QM 
((TN  *  1*0) 

<FB  *  (1.0  •  YY)) 
•  PRbF(Fa#  FS#  F) 


•  (VN  ♦  ON)  /  2*0)  •  SS  •  (VN  •  QM  •  2.0)  /  H»0 


/  (YY  •  FA) 


PRINT  85#  XL#  FA#  FB#  F#  P 
850F8RMAT  (//  15M  *‘Il*<S  LAMBDA 


•#  FtO.3  //  7H  D.F*  •#  FgtO# 


OfONrtAi  j.""--;  '  Li*  LL  -ud  .  Cn 

1AH  and#  F7.0  //  ION  F.RATie  »#  Fg.3#  5X#  3«P  •#  F7,4> 


OF  •  VN 
CC  «  TN 
08  90  1 
.  CC 
«  OF 


ON 

OF  /  2*0 
1#NF 

AL8S( 1*0 
2*0 


♦  VU)) 


CS 
OF  . 

P  »  pRBF (OF#  1000. C#  CS  /  OF) 

90  PRINT  95#  1#  XU)#  CS#  OF#  P  u.ot.uPe  .. 

950F8RMAT  (/  5H0R88T#  12#  F10.2#  lj»H  PCT.  ]jARJANj*  /f  _  -  . 

113H  CM ! -SQUARE  •#  F10.3#  5X#  6H0.F*  •#  F5,0#  5X#  3MP  »#  7*6) 

08  100  I  *  t#NV 
T(I)  ■  TU)  •  TN 
08  100  J  »  1#NQ 

100  iiirtxM  NX,  NVi  Nl. 

CALL  PRTS  (W#  NO#  NF#  5HCENT.#  Nl> 

CALL  PRTS  (C#  NV*  NF#  6WC8REL*#  Nl) 

OFH  •  TN  •  ON 

1050F8RMAT1?//  26H  UNlVARI ATE  F»TESTS*  OFB  J#  F3*0# 

16H  OFW  t#  F6*Q  ///  IBM  VARIABLE  F-RATI8#  6X#  IMP) 

08  u5  I  •  1#NV 
B  •  0*0 

08  110  J  *  l#NG 
B  *  B  ♦  S(  I# J)-*2  *  a(J) 

CC  •  T { 1 )*»2  /  TN 
F«  {(B»  CC)  *  OFH!  /  1(QU) 

P  •  pR8F(GM#  OFW#  F) 

115  PRINT  120#  I#  p#  P 
120  FORMAT  (/  16#  F12.A#  FlO.*)  . 

CALL  PRTS  (S#  NV#  N5#  6HQ  MEAN#  Ni) 

IF  (KT  *EQ.  0)  S8  T8  5 
REWIND  2  t  NT  «  TN 
08  130  1  «  i#NT 
READ  (2)  10#  (X(U)#  J  •  i#NV) 

08  125  J  *  i#NF 

125  Y{ J)  9  SCPP(X,  A*  J#  Jf  NV#  Nl) 

130  CALL  SUBS  (Y,  NF#  2H0S#  10,  N$S> 

08  T8  5  *  RETURN 

END 
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Subroutine  DMG 


This  'dispersion  matrix  generator'  subroutine  is  based  on  the  RSPACE 
program  of  Cooley  and  Lohnes  which  requires  the  input  of  the  discriminant 
function  weights,  vectors  of  group  means  and  group  dispersion  matrices,  and 
outputs  centroids  of  the  groups  and  dispersion  of  groups  in  reduced  space. 

As  our  DISCRIM  program  produces  centroids,  the  subroutine  is  necessary 
for  the  calculation  of  the  reduced  space  dispersion  matrices  only. 

The  DWTS  are  input  to  the  subroutine  in  the  call  using  the  variable  A, 
the  original  space  matrices  are  read  within  the  subroutine  from  the  rewound 
r  cch  tape  i. 

"The  r  x  r  dispersion  matrix  DDg  in  the  discriminant  space  for  group  g 
may  be  obtained  by  pre-  and  post-multiplying  the  test  space  dispersion  ma¬ 
trix  for  the  group  g  by  the  matrix  V  containing  the  discriminant  function 
(weights!  vectors  as  follows: 

DD  ,  .  =  V1.  .*  D  .  .-V,  ’’ 

gr,r)  *r,m)  g'm,m)  m,r). 


b  i 


SU3R8UYINE  DM£j(V/  NVz  Nj  <G /  NO  ^ 

DIMENSION  V(ND/N),  0(70/70  *  DMAT(70/?0)/  DDMAT(70/70) 

DIMENSION  V (NO/N ) /  0(30/30/  DMAT(30/30/DDMAT(30/30) 


PRINT  50/NV/N/KG/ND 
*>0  P8RK-AT  (10X/6I1O 
REWIND  1 
09  2  I«1/NV 

2  PRINT  7/(V{l/0/J«i*N> 

09  23  L*t*KG 

pftjlgy  Jg 

13  FOPMaT  (10X/*  DISPERSION  MATRIX  6P  INITIAL  DATA* ) 
DO  19  I*  l/NV 

19  READ  \l)  (D(1/J)/J*1/NV) 

7  P9RMAT  (10X/10P10.4) 

DO  20  I*l*N 
D9  20  J#1/NV 
DMAT(I/J>*0»0 
DO  20  K-l/NV 

20  DMAT(IzJ)  »  DMATd/J)  ♦  (V(K/1  >*D(K,J> ) 

DO  2l  I*i/N 

DO  21  J*l/N 
DDMAT(I/J)*0.0 
DO  21  KM/NV 

21  DDMATU/J)  •  DDMAT ( 1/ J)  ♦  (DMAT(I,K).V(K,J>) 

PRINT  22 

CALL  PC0S( DDMAT/ N/N/5M  DM.2/ND) 

DO  60  l*i*N 

60  PRINT  7/  (DDMATU/J)/ J*l/N)  , 

22  PORMaTUOX,*  OISPERAlON  MATRIX  OF  REDUCED  SPACE* ) 

?3  CONTINUE 
REWIND  1 

RETURN  *  END 


CLASSID  -  TPCLASS 


The  classification  routines  described  by  Cooley  and  Lohmes  in  Chapter 
seven  of  their  book  “Multivariate  Procedures  for  the  Behavioral  Sciences" 
from  the  basis  for  the  classification  of  sleep  stages  in  the  programs  CLASSID 
and  TPCLASS.  The  CLASSID  program  is  a  direct  modification  of  the  CLASSIF 
routine.  The  modified  version  is  compatible  with  the  format  of  the  training  set 
data  produced  for  the  sleep  study  and  produces  a  summary  of  the  classification 
results  for  each  group  rather  than  the  output  of  a  record  by  record  classification 
as  does  CLASSIF. 

The  program  TPCLASS  is  a  further  extension  of  this  probability  classifi¬ 
cation  system,  taking  as  data  input,  a  tape  of  unordered,  non  grouped  sleep 
measures.  The  TPCLASS  program  produces  a  plotting  of  the  probability  classi¬ 
fied  scores  versus  the  time  of  each  data  record.  It  also  provides  an  accounting 
of  the  probability  classified  score  as  compared  with  the  manually  classified 
stage  which  is  available  on  the  data  tape  used. 

In  both  programs  the  subroutine  MATINV  is  that  of  Cooley  and  Lohmes, 
a  matrix  inversion  and  determinant  calculation  by  the  Gauss-Jordan  Method 
described  in  Chapter  nine  cf  the  above  mentioned  work. 
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PRCGRAM  CUSSIDUN:»UT.8UTPUT#TaPE1/TAPE2>PU*CH) 

c  si  NO*  GROUPS#  ^VARIABLES# FIRST  GRP  NO#LAST  GRP  NO  «  512 
C  C-  NO  SUBJECTS  IN  EACH  GR8UP  ••  20P*>'0 

C  D*  D*  WEIGHTS  ••  FORMAT  1005 

C  €•  CENTR5IDS 

C  f  •  0ISPERSJ6N  MATRICES 

C  3*  FORMAT  8F  DATA  T8  DE  CLASSIFIES  —  12A6 

C  H.  Data  WITH  GR8UP  DIVIDER  CARDS  2T5  (NO  IN  GRP  ,  N8  of  GRP> 
DIMENSION  V( 50* 20)# CENT (20# 20)#  DG(50#50)#  D<20#20#20)#RATl6(2C># 
10N(20)#  A (50)#  DISC (20)#  DIF(20)#  CH! (20)#CHl$Q(20>#  Pl(20)# 
2PR8B(20)#B(50)#KLASS(20)#ATAG{20)#16(6) 

INTEGER  Fi1t(24)  #XID(1) 

DATA! I6( D# I*1#6?/10H( IX#  A5# 15#  # 10H1X#  #10H  # 

U0H(F5.3#1X)##10H*(*FA.3*)  #10HIN  *12)  / 

DIMENSION  IKAMEU2) 

1000  FORMAT  (1012) 

1001  FORMAT  (20FA.0) 

1002  FORMAT  (5E1A.7) 

1003  FORMAT  (12A6) 

1004  FORMAT  (10X#E1A«7) 

1005  FORMAT  (10X#7F10.A) 

1010  FORMaT ( IX# 12A6) 

5000  READ  1003# I NAME 

IF ( I NAME i 1 ) #EQ» INAME (2) )  STOP 

PRINT  1010# INAME 

KNT  -  0 

REWIND  1 

REWIND  2 

READ  1000#K3#M#NA#N2 
(4  •  MlNOO.J*l#M) 

ENCODE  ( 10«2# 16(3) )KG 

2  FORMAT  (I2#8H  ) 

PRINT  3# kg 

3  FORMAT  (10X#»K0«  •MlO)  „  _ 

C  VECTOR  LISTING  NO  SUBJECTS  IN  EACH  GROUP 

READ  1001#(GN(NN)#NN«1#KGJ 
PR  |  Nl^  ^0X9 

1015  FORMAT (EX#*  NUMBER  of  subjects  IN  EACH  group.) 

PRINT  1005# (GN(KN J  #NN«1#KG ) 

C  MATRIX  »F  D  WEIGHTS 
DO  12  !  *l#M 

12  READ  1009#  (V(1#J)#  J«i#N) 

C  MATRIX  0~  centroids 
DO  U  K.  1,KG 

1*  RE.  D  1005#  (CENT t l#K )#  l  •  1#N) 

£  RE>D  1005#  (CENTU#K)#K.i#KGJ 

C  DISPERSION  MATRIX 
DO  20  K*  1#KG 
DO  16  I  0  1*N 

16  READ  1005#  (OG(I#J)«  Jal#N  ) 

CALL  MATINV'*DG#N#B,  0#  DETERM  ) 

DO  18  l»l*N 
DO  18  J*1#N 
it  OU#J#K)  f  DQ( i# J) 
to  RATIO  (X)  •  GN(K)/  SORTF  (DETERM) 

PRINT  495  #(«ATI0(K»#K.1#KG) 

498  FORMAT  :ix#»  RATIO(K)  <ViOE;i2»4/lOE12*4) 

READ  1003# (FMTII )# I *1#12  ) 
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00  2Q0  UL*1*KG 

READ  202*  NSUBNG/NAGRP 

202  FORMAT  (315) 

PRINT  203/ N AGRP # NSUBNG 

203  FORMAT (lHi,»GR9UP  •  15  *  NUMBER  SUBJECTS  • 
00  303  mM«NA#N2 

303  KTAG(MH)*  MM 

PR I NT  304/ * '.IT AO ( MM  ) , MM*NA/NZ ) 

304  FORMAT { 7*t •GROUP*/ lX/ 16< ! 3/ 3X ) ) 

PRINT  305 

305  FORMAT (IX/#  ID  COUNT*) 

06  201  !»  1*16 

201  KLASSU)  •  0 

06  300  «•!* NSUBNG 

READ  FMT  ,X!D#{XU)/  I'l/M) 

KNT  a  KNT  ♦  1 


C 


DO  24  J‘1*N 
OISC(J)  •  0*0 
DO  24  In*M 

24  DISC { J )  •  0ISC(J)+  {X(I)*V(I,J)1 
DO  31  K-l/KG 
06  28  1 *i/N 

28  OIF ( j )  «  OlSC{I)*CENTtI*K) 

DO  30  J*1»N 
CHI(J)  •  0 a 0 
00  30  I«i*N 

30  CHHJ)  •  ChI(J)*(DIF(1)*DU>J*K)) 
CHISO(K)  •  0.0 

DO  3i  I “ 1* N 

31  CHISQ(K)  •  '“•USQ(K)  ♦  (DtF( t)*CH! (1 ) ) 
WRITE  (1)  anT/XI0/(CH1SQ(I)*I*1*»<6) 

33  P2  •  0.0 

DO  34  K*1#KG 

P1(K)  •  RATIO (<)  *  EXPF<-CHlSQ(K)/2.0) 

34  P2  •  P2*P1(X) 

00  36  <*1/KG 

36  PROB(K)  >  Pi (K) /P2 
TOP  «  PROB(l) 

KEEP  •  1 

08  50  I"  l/KG 

JF (Pr38 ( 1 5 .UT.T8P)G0  TO  50 


KEEP  •  ! 

TOP  a  PROS (I) 

50  continue 

KEEP  •  KEEP  ♦(NA*D 
KUASS (  KEEP)  >  KU ASS (KEEP)  *  1 
**R1NT  I6/XID/KNT/(PR03U)/I*  :G)#TPP»KEEP 
300  CONTINUE 

PRINT  800* KUASS ( NA6RP ) 

BOO  FORMAT  ( 10X* •KUASS (NAGR-,  •  •/  15) 

PCTCOP  •  ( KUASS (NAGRP) *100 *0) /NSUBNG 


PRINT  70/PCTC0R 

70  FORmaT< !OX,«PERCENT  C3RRECTL7  CUASSIF1E0  ■ 
200  CONTINUE 
GO  TO  5000 
CAUU  EXIT  6  END 


/ 
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SUBROUTINE  HAT ! NV * A, H, B, M, OETERH } 

DIMENSION  1P!V#T(50)#A(50#50)»B(50*1)*1N0EX(50*2)#PIV8T(S0) 
COMMON  PIV8T#1N0EX<1P1V9T 

EQUIVALENCE ( I ROW, JR8W ) , U C9LUM, JC8LUM ) , ( AMAX, T* SWAP ) 

CETERM  •  1.0 
15  08  20  J*1<N 
20  IP!V8T<J)«o 
30  08  550  WjN 
C  SEARCH  F9R  PIV8T  VALUE 
40  AHAX.0.0 
45  08  105  J*1«N 

so  tP( iPivert j)»i>6o#io5#40 

60  08  100  K«1,N 

70  !FUPJV8T<K>«1180ii00#74Q 

80  IF( AB9P(AMAX)*ABSPi A(Ui^) )) 85/ ICO# 100 

85  IR8W.U 

90  IC8lUm*« 

95  AMAX»A(J«K) 

100  CONTINUE 
105  CONTINUE 

110  JPIV8TCC8LUHJ.IPIV8T(IC8LUM)*1 
C  INTERCHAN3E  ROWS  T8  put  pivot  ELEMENT  on  01 ASOHAL 
130  IF ( IR8W»1C8LUH} 140#260#140 
140  0ETERM..0ETERM 
150  09  200  L«1»N 
16C  SWAP»A{ IROWiL ) 

170  A(IR9W#U»a(IC8LUM,L) 

200  A{ IC8CUH<L ) "SWAP 
205  IF(Mj260*260,210 
210  08  250  L-t#H 
220  SWAP«B(!R0W«L} 

230  B(IR8W,D»8l!C8LUK,L> 

250  B(IC8LUM#L}«8WAP 

260  iNDEXtm)  alROW 

270  I NDEX (I #  2 1 • I C8LUH 

310  PIV8T(n»A(IC8LUH|IC8LUH) 

320  DETERH«0ETERH*PlV8Tm 
C  DIVIDE  PIVOT  ROW  BY  PIVOT  ELEMENT 

330  A(!C0UUH#IC9LUH)«1«0 
340  08  350  L*1#N 

350  A ( I C9LU«j L  > • A (I C8LUM#L > /P I V8T { t ) 

355  lF(M}3a0#380#360 
360  0©  370  L*nM 

370  8 ( !C9UIM,L ) *B I!C8LUH#L ) /PIVOT ( I ) 

C  REDUCE  N8N*P!V8T  ROWS 
380  09  550  LltltN 
390  I F ( L 1 • I C8LUH ) 400 «  550 *  400 
400  T»A(Ll#IC8LUK> 

420  A(Ll»IC8LUH)«0t0 
430  09  450  L«1#N 

450  A(Ll#U)»A<LliU-A(IC8LUMiU*T 
435  IF{M)550#550#460 
460  09  500  L'l#" 

500  6{Li#U)*5(Ll#U*BnC8LUM,LH»T 
530  CONTINUE 

c  interckanqe  columns 

600  08  710  !•*  N 
610  L»N*l«I 

620  IFUNDEX(L#l)«lNOEX(L#2n630f7l0»630 


75 


630  JR8W*!nDEX<U \) 

640  JC8|»UM" ! NOEX  ( L*  2 ) 

650  08  705  <«i#N 
660  SWAP.A(K#Jr8W) 

670  AIK# jR8W>*A{</JC8UUH) 
700  A(K#gC8LUM)«SNAP 
705  CONTINUE 
710  CONTINUE 
740  RETURN 
END 
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)  R9GRAM  TPCLS( INPUT#0UTPUT#TAPe1#PUNCH) 

CCCCCC( 

cccccc 

DATA  CARDS*  A*KG#M  (312) 

8*  GNJKG)  **(20FA*0) 

C«  DWtS(V(I,J>)  (10X#6F10*A) 

D*  CENTROIDS  (CENT(I#K))  (10X#6F1Q.A> 

E*0  MAT(DG(I#J)  (lOX,7PiO.A) 

P*  NV  AL90  TITLE  (I5#6A10> 

G*  LIST  9F  VARIABLES  (2X#20I3) 

CCCCCC 

CCCCC 

C9MM9W/PP/ VM 10(300)#  NAME <  4  >  #  TPS ( 300 ) #  X ( 300 ) 

DIMENSION  T I M£S ( 6 ) # V ( 22 ) # SCORE ( 200 ) # TPSTG ( 200 ) » T AGREE ( 200 ) # 
1LVAR(22)  #XX(150)#Z(U) 

REAL  MV (33) 

ASSIGN  58  TO  NNN 
CALL  XMJT(NNN) 

CALL  8GNPLTULH.8T  #  10.0#30#50) 

1000  KKK  •  0 
PRINT  1 

1  FORMAT  (1MD 

call  class (Mv# score < 1 )  #i#kkk) 

5  CONTINUE 

1EP  «  0  •  KNT  •  0 

SUM  •  0*0 

C  NUMBER  VARIABLES  TO  BE  USED  (15)  —ANV  TITLE  *0U  WISH 
READ  2#NV#(NAM£(!)*!«1#6) 

2  FORMAT  (1X,U#0A10) 

PRINT  2#NV,(NAME(])#N1#6) 

6  FORMAT (10X#*SE0UENCE  OF  VARIABLES  TO  BE  USED*) 

READ  5/ (LVARl !)*!*1#NV) 

5  FORMAT  (2X#20|3) 

PRINT  6 

PRINT  5# (LVAR( I ) # I *1#NV) 

7Q  N  •  0 

PUNCH  300# NAME 
300  FORMAT (AAiO) 

PRINT  7 

7  F0RMaT(1HI#1X#*FILE  BEGINS  WITH  THIS  RECORD*) 

50  CONTINUE 

READ  (1>  Nt#T!MES#STAQE#aGREE#V  #D8MEGA#TAU#£ 
lF(ENOF!Le  X)  30#  11 

u  IF(n,EQ#0)  PRINT  12  #  stage#aoree#v  ,1 

12  FORMAT (IX#  F5«0#A15#2X# •PERIODS*, 

12F10.3#SF3*0/  76X#«RATES*#2F10.3#9F3.0#/1X#11F1C.A) 

N  •  N  ♦  1 
KNT  •  KNT  ♦  l 
DO  2o  L*i#NV 
KEEP  *LVAR(L) 

80  MV(L) *•  V(KEEP) 

00  21  L  •  1#U 
KK  •  L  *  9 
21  MV(KK)  •  Z(L) 

X(KNT)  •  (V(l )/1000*0) *128*0*  SUM 
XX(N)  a  X ( KNT ) 

SUM  «  X(XNT) 

SC0RE(N)*0*0 

CALL  CLASS (MV#tC0RE(N)#2#KXX) 


TPST8  (N)  •  STAGE  •  TAGREE{NJ  •  AGREE 
TPS(KHT)  •  STAGE 
VMIO(KNT)  •  SCORE (M) 

GO  TO  50 

30  IEF  •  IEP*1 
NRIGhT  •  o 

do  3i  j  •  i#n 

IF (SCORE  (<J)*EQ»TP$TG(J) )  NRIGHT  •  NRIGHT*1 
NCHK  •  10M 

:F(SCORE(W)»NC.TPSTG(J)}  NCHK  .  iOH  MISSED 

31  CONTINUE 

32  F8RMaT<1X#I3#2F10»0#2(5X,A10>#5X#£13.3) 

RIGHT  *  NRJGHT  •  COUNT  •  N 

PERCflR  •  (RIGHT/COUNT } *100.0 
print  36#jef#percor 

3*  FORMAT (iOX#*FJLE  *#12 #•  PERCENT  CORRECT  •  *#£5.2) 

PUNCH  301#IEF#PERC0R 
301  FeRMAT{lX,*F!LE**:S#Fi0.3> 

GO  T0(70#80#70#80)#|EF 
80  CONTINUE 

CALL  PICTURE(KNT) 

IF(IeF.EQ.A)  GO  TO  60 
SUM  •  0*0 

KNT  •  0  *  GO  TO  70 

58  PRINT  57,!EF,N 

57  FORMAT <1X,«XM IT  IN  CONTROL# FILE*# 15#*  RECORD*# 15) 

print  12# stage# agree# v# z 

STOP 

60  ae  to  iooo 

END 
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SUBROUTINE  CUSS  (X#PICK#  IflB#KNT) 

OIMEnSIBN  Vt 50/ 20) /CENT C20# 20}#  D&\50#50>#  &<20/20#20 '**718(30)/ 
18N(20)/  X(50)#  DI8C(20)/  0|F(20)/  CHI (20)#3{50)/Fi.ASS(2Q)#KTAQ{20) 
2#CHIS0(20)  #PR6B(20)#Pl (20) 

INTEGER  FMT<2*> 

1000  FORMAT  (312) 

1001  FBRMAT  (20F*»0) 

1002  FORMAT  (5El**7) 

1003  FBRMAT  (12A6) 

100*  FBRMAT  (10X#El*t7) 

1005  FBRMAT  (10X/7F10#*) 

IF  u00tEQ«i)  Q8  76 
10  KNT  •  KNT*  1 
06  2*  J«1#N 
0 1  SC { J )  •  0*0 
06  2A  !«1#H 
£*  DISC(J)  •  DJSC(J>* 

06  31  Kal#KQ 
06  28  I*l/N 

28  OIFU)  •  OlSC(t)*CENT(I#K) 

06  30  J»l/N 
CHI ( J)  *  0*0 
06  30  !«1/N 

30  CHI(J)  •  CHt(J)*<DIF<I  )»D(I#J,*K> ) 

CHISC(K)  c  0*0 
06  31  Ial#N 

31  CHISC(K)  •  CHI80CO  ♦  (OlF(I)*CHim) 

33  P2  •  0.0 

06  3A  K»l#XO 

ZAVE  •  EXPF(-CHI8Q(K)/2.0) 

PICK)  '•  RATIO  (X)  •  Z*VE 


100 


(X<I)*V(I#J)> 


3*  P2  •  P2  ♦Pl(K) 

06  36  K«1#X0 
36  PR60(K)  •  P1(K)/P2 
T6P  •’ PR6B{ 1 ) 

KEEP  •  1 
06  50  !•  i/KO 
l F ( FR68 ( I ) • LT » T8P ) 66  T8  50 
KEEP  «I  •  T6P  •  PR6BU) 

50  CONTINUE 

PICK  «  XEEP  »  1 
RETURN 


60  PICK  «  0.0 

PRINT '61/KN7#P2  _  _  . 

61  FORMAT  ( 10X#*£P8CH  NS.  •#!*#•  P2  •  •/0l3»5) 

RETURN 

100  KNT  •  0 

READ  1000/  *0/  H  /N1 
IF(.n6T.KQ)101/80 

101  N  *  MlNO(KQ»l#M) 

PRINT  102/KO/M#Nl 

102  FORMAT (IX# 312) 
rewind  1 

c  VECT6*  WI5TIN0  N6  SUBJECTS  |N  EACH  QR8UP 
R*.A0  1001#  (QN(NN)#NN*1#K3) 

PRINT  1005/ (QNlNN)#NN*i#KQ) 

C  MATRIX  ®F  0  WEIOHTS 
06  12  I  *l#M 

It  READ  1005/(V(1#J)/  J«1#N) 


C  MATRIX  6F  CENTRO JDS 
08  U  X*  i,KG 

1#  READ  1005*  (CENT ( J;K);  !  •  i;N) 

C  DI$P£rSJ5N  MATRIX 
OS  20  K*  l/XO 
08  16  1  ■  1#N 

X6  READ  100$;  (DG!J;J};  J»1;N> 

CALU  MAT1NV(D3;N;B;  0;  DETERM  ) 

PRINT  401#DETERM 

301  FORMAT  (20X/»DETERM»  *;E13.5) 

17  CONTINUE 

DO  18  t*i;N 
DO  is  J#1;N 

18  0(I;J;K>  •  DQ( 1/ J) 

20  RAT I 6  (K)  •  3N ( K ) /  SQRTF  (DETERM} 

PRINT  495  ;(RATI8{X)#X«tiK3) 

495  FORMAT  { IX;*  RATIO(K)  •;10El2»4/i0E12.4> 
RETURN 

80  CALL  ENDPLT  •  STOP 
END 
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disc,  -  o 

DISC,  ••  DISC,  +  X,  ■  V„  forming  t  discriminant  or  factor  score. 

DIF,  -  DISC,  -  CENT,, 

CHI.  -  0 

CHI,  -  CHI,  +  DlF,  •  D,rt  (inside  loops  over  i  and  y)  (outside  loop 

CHISQ,  -  5  o«r  i) 

CHISQ,  -  CHISQ,  +  DIF,  •  CHI,  forming the 4  group **  • 

?2  -0  * 

PI,  ~  RATIO,  ■  exponential  function  of  ( -CHISQ,/?) 

P2  -  P2  4  PI. 

PROB  -  Pl,/P2  forming  the  k  group  probability. 
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Subroutine  Picture 


Loops  1  &  2  set  up  a  matrix*  showing  how  the  computer  score  of  each 
record  compares  with  the  manual  scores  read  from  the  tape. 

VMID  is  computer  score,  TPS  is  manual  score. 

Loop  8  punches  {and  prints)  this  matrix  of  comparison  for  use  in  cor¬ 
relation  program  Correls. 

Loop  10  converts  the  computer  score  to  an  appropriate  value  for  plotting 
i.e.  stage  0  becomes  5,  lis4,  2  is  3,  3  is  2,  4  Is  1,  and  S  is  0 
Stage  5  or  REM  scores  are  converted  to  a  dotted  line  between 
stages  0  and  1 . 

The  statements  from  10  on  perform  the  plotting  according  to  the  plot 
routines  available  for  the  CDC  6600  computer  at  UT  at  Austin. 

*  This  is  a  very  compact  way  of  saving  two  scores  for  each  record  of 
data  processed  in  X(I,J)  where  X  is  count  of  the  no  of  records  manuall'- 
classified  as  I,  which  were  computer  classified  as  J. 
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SUBROUTINE  PICTURE  (NP) 

e6MM8N/PP/ VM I D  (  300 ) t NAME { 4 ) ,  TPS ( 300 ) ,  X ( 300 ) 

DIMENSION  Y(300 ) * JPT (300  5* I  SET (6/ 6) 

08  U  *  1,6 
08  1  J  •  1,6 

1  IS£T(L,J)  «  0 
*  D8  g  L  *  1,NP 

J  •  tPS(W)  ♦  J  •  K  •  VMID(U  *  1 
!F( J*E9*10)  00  TO  2 
ISET(U,K)  •  ISET(J,K)  ♦  i 

2  CONTINUE 
PRINT  11 
PRINT  76 

76  FORMAT  J15X# •MANUAU.Y  SCORED  STAGE*) 

PRINT  77 

77  FORMAT ( 10X,«  012345 

DO  8  J  *1,6 

JJ  •  J  *1 

PUNCH  99, ( I SET (K,J),K«1,6) 

99  FORMAT (613) 

8  PRINT  9#JJ,(I$ET(K,J),K  •>  1/6) 

9  FORMaT (4X/I1/5X,6(14*1a) ) 

PRINT  11 

11  FORMaTUX,//) 

DO  5  K  »  l,NP 

Y(K)  4  4.2  ft  :PT(K)  •  3 

5  CONTINUE 

DO  10  J  ■  1/NP 
VMID(J)  •  ABS(VMID(J)*5) 

IF(VmID(J).NE.O.O)  go  to  10 
VMJD(J)  •  4*0  *  IPT(J)  •  2 

10  CONTINUE 

SP  .  X(NP)/3600.0 
ENCODE ( 10*7* IXM  *X)6P 
7  F6RMaT (F5.2/5H  HRS.  ) 

ENCODE  (10/17, INPJNP 
17  FORMAT* I5,5«  PTS.  ) 

CALL  PLT(2.0,2«0/*3) 

X(NP*i)  •  0*0  •  X(NP*2)  •  5000.0 

VMIDtNP*!)  •  0#0  •  VMID(NP*2)  •  1.6 

CALL  L!NE(X,VM!0,NP,1,0#0) 

DO  90  J*  1/NP 

X(J)  •  X(J)/  X (NP  •  2) 

90  CALL  PLT(X(J),2.62,IPT(U)) 

PLACE  •  0.0  *  ICODE  *  •! 

DO  80  J  •  1*7 

CALL  SYMBOL  1 0*0* PLACE, 0*14, 9# 0*0# ICODE) 
PLACE  •  PLACE  ♦  0.625 
ICODE  •  ICODE  *1 

80  CONTINUE 

PLACE  •  0.0  •  ICODE  •  *1 
00  81  L  •  1#3B 

CALL  SYMBOL  ^LACE#  0. 0, 0. 14,9*0.0/1 COPE ) 
PLACE  •  PLACE  --  0.18 
ICODE  •  ICODE  *  1 

81  CONTINUE 

PLACE  •  *625  t  JNTEQ  •  47 
DO  9i  K  •  1*5 

CALL  SYM80L(«l*4#PLACE#Q.l05,!NTEQ,0.0#*l) 
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PLACE  «  PLACE  ♦•425  •  1NTEQ  »  1NTEQ  -  1 
9i  CBNTtNUE 

CALL  SYMBeLt-i^^i^S^OaAiliHSLEEP  STAQE#50.0>11) 

CALL  SYMB8L(t.0#-l»0#0.14#:x«AX*0.0*10> 

CALL  SYMB9Li4*0/«l*0#0«l4#!NP  ,0*0#10) 

CALL  SYMB9L(0«5#clt50#0*14iNAHe#OtO«40) 

CALL  PLT(0»OiO*0*999) 

PRINT  20*NP*VMJD{NP) 

20  PBRMAT  (10X#*PL9T  C8«PL£TED»-N8.P9!NTS»«* 18#«  LAST  f»T*«  •#£10.4} 
RETURN  9  END 


Program  Correls 


DO  60  -  sto*.  in  matrix  CODE(4,4)  the  cost  valu<>->  listed  in  XCODE(16)  of  the 

data  statement.  These  values  indicate  the  cost  of  error  of  the  computer 
classification  versus  the  manual  classification. 

100  thru  DO  3  read  in  a  classification  matrix  for  one  half  night  of  data  as  pre¬ 
pared  by  program  Tpclass. 

E>0  6  DO  7  DO  8  reduce  the  6x6  Tpclass  MATRIX  to  a  4  x  4  matrix  by  adding 
stage  1  to  stage  5REM  and  stage  3  to  stage  4. 

DO  50  Evaluate  the  e  •  re  score  matrix  by  multiplying  each  value  of  the  score 
matrix  by  the  corresponding  value  of  the  cost  matrix.  Accumulate  these 
products  in  the  variable  SUM.  If  every  decision  between  sleep  and 
awake  were  incorrect  this  sum  would  be  equal  to  the  number  of  data 
records  classified. 

DO  10-20  and  30  resort  the  4x4  computer  rs  manual  score  matrix  into  two 

vectors  of  KNT  length  where  KNT  is  no.  of  data  records  classified  -  store 
vectors  on  scratch  tape  1  to  be  read  back  for  intercorrelation  analysis. 

CALL  CORSS  (KNT,  NV)  to  perform  the  Intercorrelation  analysis.  Th;.'  routine  is 
a  modification  of  Veldmans'  CORS  subroutine  to  the  extent  that  the 
data  is  input  through  the  use  of  the  scratch  tape  1 . 

Print  50  -  Print  53  -  print  sum  as  total  error  cost.  Calculate  and  print  error  cost 
per  data  record  classified. 

DO  81  -  Print  82  -  The  total  of  the  records  on  the  diagonal  of  the  4x4  score  matrix 
represent  correct  classifcatlon.  The  percentage  of  these  correct  records 
versus  the  total  number  classified  is  calculated  ana  output. 
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MR80RAM  CORRELS  ( ! NMUT/  OUTPUT/  TAPE  I ) 

0IMENSI8N  XM(6/6)/ITlTLE(4)  /COOEU/4) 

DIMENSION  xC90E(16i 
dimension  NAME ( 8 ) 

DATA(XCOD£ ( I ) /  I  *1/ 16  I/O* / 1 •/ 1 •/ 1 •/ 1 */0*/ *25/ *5/ 1*/ *25/0*/  *23/ 
11m  *50#  *25/0*0/ 

K  •  0 

DO  60  I  •  1/4 
DO  60  v1  •  1/4 
X  •  X  ♦  l 

eODEd/J)  •  XCOOE(K) 

60  CONTINUE 

00  61  x  •  1/4 

61  MS I NT  S2/(C0DE(K/U)/L«l/4) 

52  EOrHaT (IX/ 4^5.2) 

NV  •  2 

100  READ  1/  IF IlEI /MCI/ 1  TITLE/ IF JLE2/PC2 

1  reRHAT(iX//5X«I5/riO*0/4A10/5X,I5/F10.0) 
IM(ItITLE(1)*EQ*1T!TLE(2) 1CALL  EXIT 
PRINT  2/t TITLE/ IPILEI/ IFILE2/PCI/PC2 

2  MORMaT ( 1H1/1X/4A10/*  PILE3*/2I5/2F10.3) 

REWIND  1 

DO  3  1*1/6 

READ  4/(XM(I/J)/J«l/6) 

4  FORMAT (6F3.0) 

3  MRINT  5/(XM(|/Ul/J«l/6) 

3  FORMAT (1X/6F3.0; 

DO  6  I  •  1/6 

XM{ i/2)  •  XM(I/2)4XM(l/6) 

6  XM( 1/4 )  •  XM(I,*)*XM(I/5J 
DO  7  J  *  1/4 

XMCg/J)  •  XM(2/J)  ♦  XM(6/J) 

7  XM ( 4/ J )  •  XM(4zJ>  ♦  XM(5/J) 

MRINT  70 

70  FORMAT { IX/ /// IX/ ^reduced  matrix* > 

DO  |  K  •  1/4 

a  MRlNT'5/(XM(K/U/L*l/4) 

SUM  *0*0 
00  50  I«  !/♦ 

DO  50  J»  1/A 

30  SUM  •  SUM  ♦  (XM( 1/ J)*CODEt l/Ul ) 

KNT  •  0 

DO  10  I  •  1/4 
CS  *  I 

DO  20  J  •  1/4 
SS  •  J 
L  •  XM(|,J) 

IF(U*EQ*0)QO  TO  20 
DO  90  X  *  i/L 
KNT  »  KNT  ♦  1 
WRITE  U>  XNT/CS/SS 
C  PRINT  31/XnT/CS/SS 

31  FORMAT (IX/ 13/ 2F5.0) 

30  CONTINUE 

20  CONTINUE 

10  continue 

REWIND  1 

CALL  CORSS(XNTzNV) 

PRINT  31/3UM  /KNT 
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51  P0RMaT(10X,//#10X,«ERK9R  COST  «  •F10.2i«  FOR  •  15#*  OATA  POINTS*) 

COUNT  •  KNT 
CP  •  SUM  /  COUNT 
PRINT  53#CP 

53  FORMAT U0X#*»8R  *#F10.6#*  PER  POINT*) 

DSUM  ■  0*0 
08  Sj  I  •  1#* 

81  DSUM  *  DSUM  ♦  XM(I,n 
PRINT '80# DSUM 

80  FORMAT  (10K#//#lt‘X#F6»2#*  SUM  8P  POINTS  ON  DIAGONAL*) 

PDSUM  •  ( DSUM/COUNT ) •100*0 


PRINT  82#I»DSUM 

82  F0RMaT(10:</«  OR  *#fiq.*#* 


PERCENT  ON  THE  DIAGONAL*//) 


GO  TO  100 


END 
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